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1.0  INTRODUCTION 


1.1  Impetus 

Combustion  instability  in  liquid  propellant  rocket  engines  has  historically  been  a 
major  concern  in  engine  development  programs.  Although  the  general  aspects  of  the 
combustion  process  in  liquid  propellant  rocket  combustion  chamber  are  known, 
knowledge  of  details  (e.g.,  the  rates  of  reactions  or  the  mechanisms  of  unstable  burning) 
is  incomplete.  Today  combustion  instability  can  occur  and  produce  excessive  vibration 
forces  or  heat  transfer  which  will  result  in  failure  or  corrosion  of  engine  parts, 
respectively.  The  aim  is  to  prevent  occurrence  of  this  instability  and  to  maintain  reliable 
operation.  To  meet  this  aim,  knowledge  of  the  details  governing  the  fundamental 
processes  is  required. 

1.2  Goals  and  Objectives 

Five  possible  physical  processes  associated  with  combustion  instability  are: 
primary  atomization,  secondary  atomization,  vaporization,  chemical  kinetics,  and  mixing. 
Primary  atomization  was  selected  for  the  present  study.  Strong  pressure/velocity 
disturbances  have  been  shown  to  affect  the  atomization  process  by  accelerating  the  rate  of 
jet  breakup,  thereby  changing  the  droplet  size  distribution  and  injector  spray  pattern.  The 
altered  propellant  spray  may  vaporize  at  a  different  rate  than  in  the  undisturbed  spray. 
Changes  in  droplet  size,  vaporization  rate  and  propellant  mixture  distribution  can  lead  to 
changes  in  spatial  energy  release  rate.  As  a  result,  fundamental  information  is  required 
regarding  the  breakup  of  liquid  jets  immediately  behind  a  pressure  wave. 

The  goal  of  this  study  was  to  provide  insights  associated  with  the  initial  breakup 
of  a  liquid  stream  (primary  atomization)  as  a  result  of  a  propagating  pulse.  To  meet  this 
goal,  the  objectives  of  the  study  were  to: 

1 .  Modify  an  existing  shock  tube  to  meet  the  requirements  of  the  present 
study. 

2.  Identify  the  test  liquids  and  establish  the  test  conditions  that  simulate 
engine  operating  conditions. 

3.  Conduct  the  experiments. 

4.  Analyze  the  results. 

Highly  time-resolved  images  were  employed  to  delineate  the  breakup  process. 
The  results  provide  the  understanding  needed  to  facilitate  the  development  of  improved 
atomization  models,  and  improved  steady-state  performance. 
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2.0  BACKGROUND 


2.1  Combustion  Instability  and  Its  Evolution 

Combustion  in  a  liquid  rocket  engine  is  never  perfectly  smooth;  some  fluctuations 
of  pressure,  temperature,  and  velocity  are  always  present.  When  these  fluctuations 
interact  with  the  natural  frequencies  of  the  propellant  feed  system  or  the  chamber 
acoustics,  periodic  superimposed  oscillations  occur.  In  liquid  rocket  engines  it  is  desirable 
to  minimize  the  occurrence  and  magnitude  of  chamber  pressure  oscillations  to  avoid 
catastrophic  damage  and  corrosion  of  the  combustion  chamber  due  to  large  amplitude 
pressure  waves  and  excessive  heat  transfer  rates,  respectively. 

Unstable  combustion,  combustion  instability,  displays  organized  oscillations 
occurring  at  well-defined  intervals  with  a  pressure  peak  that  may  be  maintained,  increase, 
or  die  out.  The  three  principle  types  of  combustion  instability  are  described  as  follows: 
low  frequency  {chugging}  or  system  instability  (10-200  Hz);  intermediate  frequency 
(acoustical  or  buzzing)  (200-1000  Hz);  and  high  frequency  (screaming,  or  screeching} 
(above  1000  Hz).  High  frequency  is  the  most  perplexing  and  most  common  in  the 
development  of  new  engines.  Since  energy  content  increases  with  frequency,  this  type  is 
the  most  destructive,  capable  of  rapidly  destroying  an  engine. 

High  frequency  instability  occurs  in  at  least  two  modes,  longitudinal  and 
transverse.  The  longitudinal  mode  propagates  along  axial  planes  of  the  combustion 
chamber  and  the  pressure  waves  are  reflected  at  the  injector  face  and  the  converging 
nozzle  cone.  The  transverse  modes  propagate  along  planes  perpendicular  to  the  chamber 
axis  and  can  be  broken  down  into  radial  and  tangential  modes.  Transverse  mode 
instability  predominates  in  large  liquid  rockets,  particularly  in  the  vicinity  of  the  injector. 
Thus  it  is  one  of  the  most  important  modes.  Transverse  modes  often  have  a  very 
nonlinear  characteristic  and  are  driven  predominantly  by  combustion.  Once  triggered  such 
instabilities  can  rapidly  reach  destructive  amplitudes  that  result  in  heat-transfer  rate 
increases  of  an  order  of  magnitude,  much  as  with  detonation,  causing  metal  melting  and 
wall  bum-throughs,  on  a  millisecond  time  scale. 

The  stable  combustion  field  in  a  liquid  rocket  engine  consists  of  liquid  streams 
impinging  or  shearing  due  to  high  speed  co-flowing  gas  to  form  ligaments  and  blobs  and 
finally  into  a  dense  cloud  of  droplets.  AU  of  this  liquid  bums  on  its  surface  at  a  rate 
controlled  by  its  surface  area  and  turbulent  heat  and  mass  transfer.  When  high  frequency 
instability  occurs  the  most  dramatic  change  takes  place  iu  the  poorly  atomized  and  mixed 
region  as  envisioned  by  Levine  (1965).  The  high  transverse  gas  velocity  shatters  the 
hquid,  and  particle  displacement  causes  gas-hquid  mixing,  both  resulting  in  intense  gas 
generation. 

An  understanding  of  the  origins  of  combustion  instability  could  ultimately  permit 
the  design  of  inherently  stable,  high  performance  liquid  rocket  engines  and  substantially 
reduce  the  costly  trial-and-error  development  which  is  now  necessary.  Analysis  of  data 
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could  assist  investigators  in  developing  analytical  models  which  provide  a  good 
theoretical  understanding  of  instability  mechanisms  and  in  improving  rocket  engine 
design  methods  and  procedures.  Thus,  current  instrumentation  and  diagnostic  capabilities 
enable  a  new  approach  to  the  control  of  combustion  instability;  investigation  and  eventual 
control  of  the  basic  combustion  instability  mechanisms. 


2.2  Prevention  of  Combustion  Instability:  The  Motivation 

Despite  advances  in  understanding  and  modeling  combustion  instability 
processes,  conducting  stability  rating  programs,  and  designing  injectors  and  combustion 
systems  to  avoid  instability,  unstable  operation  is  still  encountered.  Combustion 
instability  thus  continues  to  be  a  concern,  particularly  now  that  a  new  round  of  liquid 
rocket  engine  development  is  imminent. 

Current  approaches  for  passive  control  of  combustion  instability  rely  on  damping 
out  pressure  oscillations  by  the  use  of  baffles,  resonators,  or  acoustic  liners.  Active 
control  of  combustion  instabilities  involves  the  control  of  energy  addition/extraction  of 
combustion/flow  processes.  Advances  in  fast-response,  nonintrusive  instrumentation, 
measurement,  and  diagnostic  techniques  now  make  possible  the  investigation  of  many 
combustion  processes  which  could  participate  in  liquid  rocket  engine  instability,  and 
thereby  (1)  create  the  possibility  of  tmderstanding  the  origins  of  combustion  instability, 
and  (2)  permit  ultimately  the  design  of  inherently  stable,  high  performance  liquid  rocket 
engines. 

2.3  Candidates  for  Instability  Mechanisms 

Continual  effort  over  the  past  several  decades  has  identified  numerous  candidates 
for  instability  mechanisms.  A  JANNAF  sponsored  workshop  on  Liquid  Rocket  Engine 
Combustion  Driven  Instability  Mechanisms  resulted  in  identification  of  key  mechanisms 
and  prioritization  of  research  needs  associated  with  these  mechanisms.  Five  categories  of 
basic  combustion  physics  associated  with  combustion  instabilities  were  identified 
(Jensen,  1989); 


•  Primary  atomization 

•  Secondary  atomization 

•  Vaporization 

•  Chemical  kinetics 

•  Mixing 

Of  these  five,  the  category  assessed  to  have  a  great  potential  for  participation  in  an 
instability,  and  therefore  selected  for  investigation,  is  the  effect  of  transverse  mode 
instability  on  primary  atomization  of  the  liquid  jet.  The  committee  concluded  that 
experimental  studies  of  both  atomization  and  spray  characteristics  in  steady  and  simulated 
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instability  conditions  should  be  emphasized.  Complementary  cold-flow  simulations  along 
with  limited  combustion  testing  should  be  considered  to  facihtate  extrapolation  of  models 
based  on  more  extensive  cold  flow  data.  Objectives  are  detailed  atomization  rates  and 
drop  sizes  and  gas  phase  parameters  to  anchor  mechanism  models.  Sub  and  supercritical 
conditions  were  also  recommended  for  consideration.  The  subsequent  discussion  explains 
the  selection  of  primary  atomization. 

Three  fundamental  multi-phase  conditions  exist  within  the  combustion  chamber: 
1.)  liquid  core/ligament  regime,  2.)  dense  spray  regime,  and  3.)  dilute  spray  regime.  There 
exists  multiple-processes  within  each  regime.  Dominant  processes  vary  from  regime  to 
regime  and  form  distinct  zones  within  each  regime.  The  transient  nature  of  the  dispersed 
interfacial  boundary  between  dense  and  dilute  spray  regimes  is  a  critical  factor  in  terms  of 
flame-stabilization. 

Each  process  has  a  characteristic  time.  Variation  in  pressure,  temperature,  and 
velocity  will  affect  the  atomization  process  and  the  characteristic  times  and  rates 
associated  with  those  processes.  If  one  of  the  characteristic  times  is  substantially  longer 
than  the  others,  that  particular  process  can  be  rate-controlhng.  This  can  introduce 
dynamic  or  frequency-dependent  behavior  which  will  amplify  the  pressure  and  velocity 
waves  to  considerable  strength.  Since  these  pressure/velocity  fluctuations  have  been 
shown  to  accelerate  the  breakup  mechanism,  the  atomization  process  wUl  tend  to  be  a 
dynamic  feedback  behavior. 

2.4  Primary  Atomization 

Primary  atomization  in  a  liquid  rocket  engine  arises  from  the  breakup  of  liquid  jets 
or  sheets,  depending  on  the  type  of  injection  element  employed.  In  this  atomization 
mechanism  initial  separation  of  large  drops  and  ligaments  from  the  main  body  of  the  hquid 
occurs  under  very  high  inertial  forces.  At  the  high  injection  velocities  used,  Rayleigh 
breakup  is  not  relevant,  and  the  Kelvin-Hehnholtz  liquid  instability  is  of  greater  interest. 
In  this  process,  surface  waves  form  and  grow,  causing  distortion  and  eventual  breakup  of 
the  fluid  surface. 

2.5  Jet  Breakup  and  Penetration 

In  the  breakup  of  a  liquid  jet,  the  jet  is  suddenly  subjected  to  large  aerodynamic 
forces  that  result  in  jet  breakup.  Investigations  of  liquid  jet  breakup  by  a  shock  wave  have 
demonstrated  that  the  jet  is  broken  by  the  high  velocity  gas  flow  behind  the  shock  front 
(Morrell,  1963).  The  flow  behind  the  shock  wave,  moving  transverse  to  the  liquid  jet, 
causes  disintegration  to  occur  simultaneously  at  all  points  along  its  length. 

In  the  present  study,  the  initial  response  of  the  liquid  to  the  high  speed  gas  flow 
occurs  due  to  the  external  pressure  distribution  around  the  cylindrical  liquid  jet.  The  flow 
around  the  jet  acts  like  flow  around  a  solid  cylinder.  For  example.  Figure  1  illustrates  the 
pressure  coefficient  distribution  around  a  solid  cylinder.  A  similar  pressure  distribution 


4 


creates  an  adverse  pressure  gradient  within  the  liquid,  resulting  in  flattening  of  the  liquid 
jet  transverse  to  the  flow  of  gas.  At  the  outer  edges  of  this  flattened  sheet,  drops  and 
ligaments  are  tom  off  when  the  tangential  shear  stress  overcomes  the  surface  tension  of 
the  fluid.  The  rate  of  breakup  is  assumed  to  be  controlled  by  the  distortion  rate  of  the 
liquid  cross  section.  The  model  in  Figure  2  illustrates  the  proposed  primary  mechanism  of 
breakup  of  a  liquid  jet  in  a  cross  flow  of  gas.  High  viscosity  may  retard  the  distortion,  but 
the  controlling  step  should  remain  the  same.  In  the  present  study,  the  time  required  for 
complete  breakup  to  occur  was  judged  visually  from  stUl  frames  and  high  speed 
photographs. 

Experimental  observations  of  jet  behavior  after  shock  passage  have  included  the 
time  tj  for  the  jet  to  shatter,  its  drag-induced  displacement  x  from  the  initial  location,  the 
height  h  of  jet  penetration,  and  the  coordinates  x,  y  of  the  jet  (Reinecke,  1977;  Kush,  and 
Schetz,  1972;  Reichenback,  and  Horn,  1971).  These  quantities  generally  have  been 
nondimensionalized  by  the  diameter  Dj  and  density  py  of  the  jet  and  by  the  air  density  pg 
and  speed  Vg  behind  the  shock  wave  resulting  in  the  dimensionless  quantity, 

X  =  xlDj  (2-1) 


Figure  1.  Pressure  Coefficient  Distribution  Around  a  Solid  Cylinder 
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Figure  2.  Primary  Mechanism  of  Liquid  Jet  Breakup 
The  penetration  commonly  is  nondimensionalized  and  correlates  as  (Reinecke, 

1977): 


(2-2) 


where  D  is  the  jet  initial  diameter,  M  and  p  are  the  freestream  Mach  number  and  static 
pressure,  and  p^j  is  the  jet  plenum  pressure.  Since  the  liquid  jet  is  usually  underexpanded, 
Eq.  (2-2)  can  also  be  written  as: 


h/Dj  =  K 


ypj^j 


(2-3) 


where  pj  Vj  and  are  the  jet  and  airstream  dynamic  pressure,  respectively,  and  y  is 

the  specific  heat  ratio  of  the  air.  Typical  values  for  K  are  about  7.0  in  experiments  where 
the  Weber  number  was  generally  above  500. 

Assuming  that  the  jet  element  moves  normal  to  the  air  velocity  at  constant  speed 
Vj,  the  distance  h  that  the  jet  penetrates  into  the  flow  is  the  product  of  the  time  tj  required 
to  shatter  the  jet  element  times  the  speed  of  the  jet,  that  is  h=tjVj.  From  this  relation  and 
from  Equation  (2-3),  the  dimensionless  time  for  breakup  is  derived: 
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(2-4) 


(forY=1.4,K=7.0) 


The  downstream  displacement  of  the  jet  can  also  be  calculated.  Values  of  the 
downstream  displacement  of  the  jet  element  as  a  function  of  its  time  after  injection  y/Vj, 
where  x  and  y  are  the  coordinates  of  the  jet,  were  determined  from  each  frame  captured 
at  33  ^is  intervals.  A  plot  of  these  data  in  the  nondimensional  form  of  X  =  x/Dj  as  a 


tV  p 

function  of  T  =  — -  —  was  generated  to  compare 


P; 


agreement  of  the  empirical  coefficients. 


the  results  for  correlation  and 


2.6  Summary  of  Relevant  Research 

Combustion  instabihty  in  hquid  propellant  rocket  engines  has  historically  been  a 
major  concern  in  engine  development  programs  and  was  the  subject  of  extensive  research 
during  the  60's  and  early  70’s  (Hanje  and  Reardon,  1972).  Despite  advances  in  modehng, 
stability  rating  techniques,  stability  aids,  and  injector  design,  unstable  operation  is  still 
encountered.  Combustion  instability  continues  to  be  a  concern,  particularly  now  that  a 
new  round  of  liquid  rocket  engine  development  is  imminent. 

Current  approaches  for  passive  control  of  combustion  instabihty  rely  on  stabhity 
aids  such  as  baffles  and  acoustic  cavities,  and  injector  design,  although  active  control  of 
combustion  instabhities  has  also  been  recently  proposed.  This  “band-aid”  approach  is 
necessitated  by  a  lack  of  fundamental  imderstanding  of  how  combustion  instabhities  arise, 
grow,  and  sustain.  An  understanding  of  the  origins  of  combustion  instabihty  could 
ultimately  permit  the  design  of  inherently  stable,  high  performance  hquid  rocket  engines 
without  much  of  the  costly  trial-and-error  development  which  is  now  necessary. 

Experimental  analysis  of  basic  instabihty  mechanisms  could  assist  investigators  in 
developing  analytical  models  to  provide  a  good  theoretical  imderstanding  of  how  these 
mechanisms  cause  instabihty  and  ultimately  provide  a  design  tool  for  improving  rocket 
engine  design.  Recent  advances  in  fast-response,  non-intrusive  diagnostic  techniques  now 
make  possible  the  investigation  of  many  combustion  processes  which  could  participate  in 
hquid  rocket  engine  combustion  instability. 

Over  the  past  several  decades  candidates  for  instabihty  mechanisms  have  been 
identified  by  researchers  (primary  atomization,  secondary  atomization,  vaporization, 
chemical  kinetics,  and  mixing).  Of  these  five,  secondary  atomization  is  the  subject  of  one 
Air  Force  program  (Cox,  1990).  The  current  effort,  also  Air  Force  sponsored,  is  designed 
to  assess  the  effect  of  transverse  mode  instabihty  on  primary  atomization  of  the  hquid  jet 
common  to  most  rocket  engine  injection  systems.  The  aforementioned  JANNAF 
sponsored  committee  concluded  that  experimental  studies  of  both  atomization  and  spray 
characteristics  in  steady  and  simulated  instabihty  conditions  should  be  emphasized  and  that 
complementary  cold-flow  simulations  along  with  limited  combustion  testing  should  be 
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considered  to  facilitate  extrapolation  of  models  based  on  more  extensive  cold  flow  data. 
Targeted  data  included  detailed  atomization  rates  and  drop  sizes  and  gas  phase  parameters 
to  anchor  mechanism  models. 

In  a  liquid  rocket  engine,  primary  atomization  is  obtained  by  impinging  liquid 
streams  or  shearing  a  liquid  stream  with  a  high  speed  co-flowing  gas.  In  a  transverse 
spinning  mode  instability,  wave  induced  primary  atomization  of  the  propellant  streams 
causes  greatly  increased  generation  of  droplets  and  thus  surface  area  available  for 
vaporization.  The  convective  flow  associated  with  the  wave  passing  enhances  the 
evaporation  rate,  gas  phase  mixing,  and  reaction  rates.  Additionally,  reduced  penetration 
of  the  spray  into  the  combustion  chamber  occurs,  increasing  the  concentration  of  reactive 
components  near  the  injector  and  increasing  the  energy  release  rate  in  that  vicinity.  If  the 
expanding  gases  created  during  energy  release  increase  the  pressure  amplitude  of  the 
passing  wave  or  the  amplitude  of  an  ensuing  wave  and  overcome  damping  processes,  an 
mstability  can  occur. 

When  high  frequency  instability  occurs,  the  most  dramatic  change  probably  takes 
place  in  the  poorly  atomized  and  mixed  region.  As  envisioned  by  Levine  (1965),  the  high 
transverse  gas  velocity  shatters  the  liquid,  and  particle  displacement  causes  gas-liquid 
mixing,  both  resulting  in  intense  gas  generation.  At  low  amplitudes  these  effects  may 
damp  as  well  as  drive,  but  at  high  amplitude  the  net  effect  occurs  at  the  proper  phase 
relationship  to  sustain  the  high  frequency  instability. 

Considerable  attention  has  been  given  to  the  physical  process  of  injection  and  its 
role  in  the  control  of  the  rate  and/or  completeness  of  combustion  [Clark  (1964);  Heideman 
(1968);  Buffam  (1967);  Huynh;  (1992);  Sitkei  (1963)].  Most  of  the  interest  has  been 
concentrated  on  the  drop  sizes  produced  by  different  types  of  injectors  under  various 
environmental  conditions.  Mass  distribution  and  drop  size  measurements  have  been  made 
for  the  breakup  of  fuel  jets  in  a  crosscurrent  flow  of  gas  (Clark,  1964).  Others  have 
studied  the  behavior  of  drops  and  sprays  in  the  gas  flow  in  a  shock  tube  or  blast  gas  or  in 
different  types  of  tunnels.  An  analysis  of  critical  conditions  for  jet  breakup  indicates  that 
there  is  a  predictable  velocity  at  which  breakup  commences  and  that  this  velocity  is  a 
function  of  the  liquid  dimensions  and  properties  and  the  duration  of  the  flow  (Morrell, 
1961).  It  appears  that  jet  shattering  may  be  a  rate  process  that  can  account  for  the  growth 
and  propagation  of  a  wave  in  a  liquid-fuel  combustion  system  that  is  vaporization  limited. 
In  oscillatory  combustion,  Clark  (1960)  identified  a  possible  mechanism  for  amplifying  the 
instability  through  the  increase  in  energy  release  rate  caused  by  more  rapid  atomization  in 
the  wake  of  the  pressure  wave. 

Missing  are  high-resolution,  high-speed  photography  to  delineate  the  processes 
associated  with  the  breakup  of  a  liquid  stream  by  a  transverse  shock  wave.  Previous 
images  of  jet  breakup  in  a  crossflow  have  been  limited  to  low  time  resolutions  on  the  order 
of  milliseconds.  However,  the  different  regimes  the  jet  undergoes  before  complete 
breakup  occur  much  before  1  millisecond  in  high  speed  gas  flows.  In  the  present  study,  a 
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cylindrical  liquid  jet  formed  by  a  simple  orifice  acted  upon  by  a  transverse  shock  wave  was 
investigated  experimentally  by  capturing  images  at  a  33  |xs  time  interval. 

In  past  studies,  the  liquid  jet  was  subjected  to  a  constant  gas  flow  velocity.  Based 
on  the  literature  review,  no  studies  have  been  conducted  to  study  the  behavior  of  a  hquid 
jet  in  an  exponentially  decaying  velocity  field  (N-wave).  The  behavior  of  drops  have  been 
studied  in  an  accelerating  and  decelerating  flowfield,  and  it  was  shown  that  drops  do 
respond  to  the  deceleration  forces.  It  was  with  this  intent  that  the  behavior  of  hquid  jet  in 
such  an  environment  was  studied. 

In  addition,  past  experiments  which  utihzed  photographic  images  to  delineate  the 
atomization  were  limited  to  large  hquid  to  gas  momentum  ratios  (q  >  1).  And  it  was 
shown  that  the  jet  first  deforms  before  any  breakup  occurs.  However,  this  behavior  was 
only  studied  at  high  momentum  ratios.  In  this  study,  the  simulation  of  a  hquid  rocket 
engine  operating  condition  was  of  paramount  concern.  In  order  to  do  so,  low  value  of 
momentum  ratios  had  to  be  obtained  (q  «  1)  to  accurately  delineate  the  breakup  process. 

Insight  into  this  mechanism  is  sought  to  assess  its  potential  for  causing  combustion 
instabihty  in  a  hquid  rocket  engine.  Atomization  and  breakup  rate  for  any  point  in  the  jet 
can  be  calculated  from  the  experimental  data. 


3.0  APPROACH 

Task  1.  Preparation  of  Hardware  to  Simulate  Transverse  Wave 

This  task  addresses  the  modification  of  an  existing  shock  tube  to  allow  the 
presentation  of  a  hquid  jet.  The  shock  tube  should  repeatably  produce  a  pressure  pulse  of 
the  correct  amphtude,  duration,  and  character.  In  addition,  the  shock  tube  should  include: 
(1)  a  pressurized  injector  housing  to  dehver  the  hquid  simulant  into  the  bore  of  the  shock 
tube  at  both  atmospheric  and  high  pressure  conditions,  (2)  a  fluid  dehvery  system 
consisting  of  a  pulsation-  and  vibration-free  pressure  feed  system  for  atmospheric  test 
section  conditions,  and  a  gear-pump  feed  system  for  high  pressure  experiments,  and  (3) 
stream  injectors  to  dehver  a  laminar  hquid  jet  column  into  the  shock  tube.  Finahy  a 
selection  of  proper  hquid  simulants  of  hquid  oxygen  (LOX)  is  required  to  establish  of  test 
conditions  that  simulate  rocket  engine  jet  breakup.  Surface  tension  and  viscosity 
properties  of  the  simulants  need  to  be  matched  with  LOX  properties  in  order  to  simulate 
the  breakup  phenomena  as  well  as  possible. 

Task  2.  Development  of  Diagnostics  and  Data  Acquisition  System 

In  order  to  investigate  the  shock-induced  stream  breakup,  a  fast-response  imaging 
and  diagnostic  technique  must  be  developed.  High  speed  shadowgraph  images  (30,000 
frames/sec)  of  the  transient  behavior  of  the  hquid  jet  are  required  using  a  high  intensity, 
short  duration,  back-hght  source.  Pressure  transducers  are  needed  upstream  of  the  hquid 
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jet  to  sense  the  approaching  shock  and  initiate  the  cinematography  sequence  via  a 
microprocessor  controller. 

A  data  acquisition  system  is  also  required  to  acquire,  display,  analyze,  and  store 
data  from  the  pressure  transducers  signals.  The  resulting  pressure  plot  histogram  is 
needed  to  calculate  many  of  the  gas  phase  properties  (i.e.,  velocity,  density,  temperature, 
and  viscosity).  This  information  can  then  be  used  to  determine  instantaneous  values  of 
nondimensional  parameters  that  are  deemed  important  in  the  breakup  process. 

Task  3.  Development  and  Implementation  of  a  Test  Matrix 

Experimental  modeUng  of  the  atomization  process  inside  a  high  pressure  hquid 
rocket  combustion  chamber  requires  the  matching  of  certain  fluid  properties  and 
dimensionless  parameters  deemed  important  in  this  atomization  process  of  the  liquid 
oxygen.  A  test  matrix  is  needed  to  identify  initial  shock  tube  and  stream  injector 
conditions,  for  a  particular  hquid  simulant,  that  simultaneously  matches  most  of  the 
atomization  parameters  in  a  particular  hquid  rocket  engine.  The  engine  chosen  to  be 
studied  is  the  LOX/H2  Space  Shuttle  Main  Engine  (SSME). 

Task  4.  Data  Analysis 

The  breakup  mechanism  and  wave  interaction  must  be  analyzed  and  quahtative  and 
quantitative  data  analysis  employed  in  the  investigation  of  hquid  jet  breakup.  Captured 
images,  for  example,  can  provide  information  on  distortion  of  the  hquid  round  jet,  the 
formation  of  axial  surface  waves,  jet  displacement,  drop  size,  etc.  Quahtative  analysis 
consists  of  studying  this  general  behavior  of  the  jet  at  various  operating  conditions. 
Quantitative  analysis  consists  of  measurements  extracted  directly  from  the  photographs 
and/or  from  digitahy  scanned  images. 


4.0  EXPERIMENT 


4.1  Overview 

To  investigate  the  stream  breakup  a  planar  pressure  wave  was  passed  traversely 
through  a  round  hquid  jet.  A  shock  tube  was  the  vehicle  used  to  simulate  the  passage  of  a 
transverse  instabhity  wave.  The  plane  shock  wave  was  produced  by  the  sudden  bursting 
of  a  diaphragm  which  separates  a  gas  at  high  pressure  from  one  at  lower  pressures.  After 
the  bursting  of  the  diaphragm,  a  compression  wave  forms  in  the  low-pressure  gas,  this 
rapidly  steepening  to  form  a  shock  front.  At  the  same  time  a  series  of  expansion  waves 
propagate  in  the  high-pressure  gas.  Curves  of  p  versus  x  and  T  versus  x  shortly  after 
rupture  of  the  diaphragm  are  shown  in  Figure  3  (John,  1984). 

The  bursting  of  the  diaphragm  generates  a  shock  wave,  which  propagates  into  the 
expansion  chamber,  and  an  expansion  wave,  which  propagates  into  the  compression 
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chamber,  as  shown  in  Figure  3.  The  condition  of  the  fluid  which  is  traversed  by  the  shock 
is  denoted  by  (2),  and  that  of  the  fluid  traversed  by  the  expansion  wave  is  denoted  by  (3). 
The  interface  between  regions  2  and  3  is  called  the  contact  surface.  It  marks  the  boundary 
between  the  fluids  which  were  initially  on  either  side  of  the  diaphragm.  On  either  side  of 
the  contact  surface,  the  temperatures,  T2  and  T3,  may  be  different  (only  if  conduction  is 
neglected),  but  it  is  necessary  that  the  pressure  and  fluid  velocity  be  the  same,  that  is,  p2  = 
P3.  Therefore,  the  contact  surface  travels  with  the  speed  of  the  gas  flow  behind  the  shock. 
During  the  time  that  the  images  were  captured  (1  msec  after  the  passage  of  the  shock)  the 
contact  surface  did  not  have  enough  time  to  reach  the  test  section  for  the  gas  density 
variation  to  be  an  issue.  In  addition,  the  driven  section  following  the  test  section  was 
suf^cigntly  long  to  not  allow  the  reflected  shock  to  reach  the  test  section  during  the 
cinematography. 

A  double  diaphragm  high-pressure  shock  tube  was  used  to  provide  the  pulse 
required  to  perturb  a  Uquid  jet  and  induce  the  primary  atomization  instability  mechanism 
(Figure  4).  The  shock  tube  is  square  in  cross-section  (2"  x  2")  and  capable  of  holding 
pressures  to  1000  psia  in  the  driven  section  and  3000  psia  in  the  driver  section.  The 
double  diaphragm  technique  for  pulse  generation  can  repeatedly  provide  on-demand, 
individual  pulses  of  the  correct  character,  duration,  and  amplitude. 

Three  different  driver  section  configurations  were  used.  For  generating  step  waves 
a  36”  cylindrical  driver  was  used.  For  generating  N-waves,  a  conical  driver  section  was 
used. 


The  technique  used  to  produce  N-waves  is  based  on  the  spherical-balloon  problem 
(Temkin,  1982).  Here  a  pressurized,  spherical  region  is  allowed  to  expand  suddenly.  As 
was  found  (Lamb,  1925),  the  expansion  produces  an  N-wave.  By  analogy,  a  long  cone 
with  a  pressurized  region  near  the  vertex  can  also  be  used  to  produce  such  waves.  For 
this  reason  a  conical  driver  section  was  selected,  while  the  driven  section  was  of  constant 
cross-sectional  area.  A  36”  and  an  11”  conical  driver  section  were  used  to  produce  two 
different  gas  flow  duration.  The  36”  conical  driver  produced  a  wave  5.5  msec  in  duration, 
whereas  the  11”  driver  produced  a  2.5  msec  wave  duration.  These  times  closely 
correspond  to  the  acoustic  limit  t^  =  2X(^a  where  Xq  is  the  length  of  the  pressurized  cone, 
and  a  is  the  speed  of  sound  in  the  low-pressure  region. 

Figures  5a  and  5b  show  a  typical  pressure-amplitude- versus-time  record  for  a  step 
wave  and  an  N-wave,  respectively,  obtained  with  the  facility. 
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Figure  3.  Pressure  and  Temperature  as  a  Function  of  x 


4.2  Stream  Injector 

The  liquid  is  injected  into  the  bore  of  the  shock  tube  via  a  pressure  feed  or  a  pump 
system.  This  system  of  fluid  delivery  is  pulsation  and  vibration  free.  A  hypodermic  needle 
and  precision-hole  stainless  steel  tubing  were  used  as  an  injector  orifice.  Higher  pressure 
testing  requires  the  use  of  a  gear  pump  to  feed  the  liquid  into  the  shock  tube.  The  injector 
is  placed  inside  a  high  pressure  housing  with  an  internal  traverse  system  to  vary  the 
penetration  length  viewed  from  the  side  optical  access. 


4.3  Liquid  Simulant 

Most  important  to  the  establishment  of  the  subcritical  test  conditions  is  the 
selection  of  the  proper  liquid  simulants  (Table  1).  Since  viscosity  and  surface  tension  need 
to  be  matched  with  the  LOX  properties,  these  are  the  variables  which  require  the  most 
scrutiny.  Methanol  was  used  for  preliminary  chamber  checkout  tests.  Freon-22  was 
chosen  as  the  simulant  of  choice  for  subcritical  tests.  For  supercritical  testing  Freon  116 
and  Freon  503  have  been  identified  to  be  the  simulants  of  choice  due  to  their  low  critical 
temperature  and  pressure  (Table  2). 
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Liquid  from  facility 


Figure  5.  Typical  Pressure- versus-Time  profile  (a)  Step  wave,  (b)  N-wave 


Table  1.  Subcritical  Simulant  Liquids 


Liquid 

a  O^g/s^) 

Freon  -  22  @  250C 

1194 

0.008 

1.98 

Methanol  @  250C  latm 

791 

0.0226 

5.62 

Water  @  5^C  1  atm 

1000 

0.0731 

9.67 

LOX  @  170OR 

1153 

0.0121 

2.20 

Table  2.  Supercritical  Simulant  Liquids 


Fluid 

Ter  TO 

Prr(PSIA) 

■99^91 

a  -  (kg/s2) 

Freon  116 

68 

432 

1587  @-730C 

0.016  @-730C 

3.00  @  -25OC 

Freon  503 

67 

632 

1223  @  -30OC 

0.0061  @  -30OC 

1.44  @  -30OC 

LOX  @170OR 

-183 

732 

1153 

0.0121 

2.20 
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4.4  Facility 

A  schematic  of  the  flow  facility  is  shown  in  Figure  6.  A  control  panel  meters  the 
flow  into  the  shock  tube  via  series  of  regulators  and  solenoid  valves.  The  firing  sequence 
of  the  shock  wave  is  as  follows.  Solenoid  valves  SI,  S2,  S3  are  opened  and  S4  closed. 
Volumes  VI,  V2,  and  V3  are  filled  to  the  desired  pressure  of  V3.  S3  is  closed  and  VI, 
and  V2  are  filled  up  to  the  desired  pressure  of  V2.  S2  is  then  closed  and  VI  is  filled  to  the 
desired  pressure  in  the  driver  section  and  then  SI  is  closed.  To  initiate  the  shock  S4 
solenoid  is  opened  to  exhaust  the  gas  in  the  intermediate  section  (V2)  into  the  driven 
section  (V3).  This  causes  a  high  enough  pressure  differential  across  the  first  and  the 
subsequent  diaphragm  to  cause  bursting.  This  method  of  shock  initiation  produces  an  on- 
demand  repeatable  shock  wave  of  correct  character,  duration,  and  amplitude. 

4.5  Diagnostics 

In  order  to  investigate  the  shock-induced  stream  breakup,  a  fast-response  imaging 
and  diagnostic  technique  has  been  utilized  to  investigate  the  fluid  mechanics  associated 
with  the  combustion  processes  which  could  participate  in  the  evolution  and  development 
of  combustion  instability.  Figure  7  illustrates  the  diagnostic  set-up.  A  planar  pressure 
wave  passes  transversely  through  a  Uquid  round  jet.  Two  high-frequency  piezoelectric 
pressure  transducers  are  positioned  on  either  side  of  the  injection  port  on  the  test  section. 
These  miniature  dynamic  transducers  are  specifically  designed  for  shock  tube 
measurements  where  very  high  frequency  response  is  required.  A  digital  delay  /  pulse 
generator  (Stanford  Research  Systems  Model  DG535)  was  used  to  issue  burst  signals  to 
the  laser  for  still  image  photography  and  high  speed  cinematography. 

Hieh  Speed  Camera 


Screening  tests  were  conducted  using  a  35  mm  Pentax  camera  to  capture  stiU 
images  of  jet  breakup  at  different  time  intervals.  Each  image  in  the  breakup  sequence  is  a 
result  of  different  tests  conducted  under  the  same  shock  tube  condition.  Results  of 
pressure  plots  indicate  very  good  shock  repeatability  to  within  0. 1  psi.  High  speed  images 
of  the  transient  breakup  were  captured  with  a  Cordin  (Model  321)  rotating  drum  camera. 
The  drum  camera  operates  by  the  principle  of  a  rotating  drum  with  aim  length  35  mm 
film  on  the  track.  A  camera  controller  is  utilized  to  supply  the  proper  current  load  to  the 
drum  motor.  The  drum  can  rotate  at  a  maximum  rate  of  300  revolutions/sec  (rps).  As  tiie 
film  rotates  with  the  drum  the  shutter  is  opened  when  the  shock  is  initiated,  thus  allowing 
pulsating  light  to  be  directed  through  the  camera  optics  onto  the  film.  Each  pulse  of  light 
will  direct  a  discrete  image  on  the  film.  An  Infimty  Model  K2  long  distance  microscope 
was  used  to  resolve  the  breakup  phenomena  accurately. 
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CAMERA 

CONTROLLER 


Copper  Vapor  Laser 


The  Copper  Vapor  Laser  (CVL)  high  intensity  visible  light  source  offers  many 
benefits  in  high  speed  liquid  breakup  visualization.  The  ultra  short  pulse  width  (20  ns) 
freezes  the  wave  induced  motion  of  the  liquid  jet.  In  the  screening  experiments  the  CVL 
was  in  the  ‘single’  mode  operation  for  still  image  photography.  In  high  speed 
cinematography  the  "Multishot"  operation  was  selected  (30  kHz)  and  was  synchronized  to 
the  Cordin  camera  via  an  'n'-shot  controller.  The  controller  is  a  microprocessor  based 
device  which  allows  the  laser  to  be  operated  in  burst  mode.  The  controller  operates  the 
laser  at  a  free  run  frequency.  The  CVL  is  then  switched  off  momentarily  at  which  point  a 
burst  of  'n'  pulses  at  a  set  burst  frequency  are  fired.  After  a  second  time  delay  the 
controller  continues  to  operate  the  laser  at  the  free  run  frequency.  A  shadowgraph 
technique  was  implemented  for  capturing  the  images.  The  camera  was  positioned  in  front 
of  the  side  viewport  optic  with  the  laser  penetrating  through  the  two  test  section 
viewports  backlighting  the  liquid  jet.  The  timing  sequence  is  illustrated  in  Figure  8.  The 
trigger  signal  from  the  shock-initiating  solenoid  valve  initiates  a  sequence  of  events.  The 
solenoid  signal  is  sent  to  the  ‘n’  shot  controller  where  a  signal  is  issued  to  the  camera 
shutter,  allowing  it  to  open.  At  this  time,  a  signal  is  also  issued  to  the  laser,  turning  it  off 
momentarily,  to  prevent  film  exposure  during  the  opening  of  the  shutter.  After  the  shock 
wave  has  propagated  downstream,  it  triggers  the  first  pressure  transducer.  The  signal  is 
sent  to  a  delay  generator  where  after  a  set  delay,  a  TTL  signal  is  sent  to  the  laser  .signaling 
it  to  issue  a  burst  of  ‘n’  pulses.  Each  pulse  places  a  discrete  image  onto  the  image  plane. 
After  the  bursts  are  complete,  the  shutter  is  signaled  to  close.  And  after  a  set  time  delay, 
the  laser  reverts  back  to  its  normal  free  run  frequency. 

A  Kodak  EktaPro  EM  Motion  Analyzer  System  (Model  1012)  was  also  used  with 
room  light  to  capmre  high  speed  video  images  of  the  jet  breakup  as  viewed  from  the  side 
and  end  viewport. 

4.6  Data  Acquisition  System 

The  TestLab  multi-channel  waveform  analyzer  (Model  2500)  was  used  to  acquire, 
display,  analyze,  and  store  data  from  the  pressure  transducers  and  the  burst  mode  of  the 
CVL.  The  TestLab  provides  spreadsheet-style  management  and  set  ups  of  data  files  which 
are  transferable  to  an  IBM  PC.  The  transient  pressure  plots  are  plotted  and  are  used  to 
calculate  the  time  variant  gas  flow  velocity,  Weg  &  Reg,  liquid/gas  dynamic  pressure 
ratio,  and  viscosity  ratio.  A  photodiode  is  utilized  to  acquire  signals  from  the  laser  pulses 
to  coordinate  the  time  each  image  was  captured  with  its  instantaneous  pressure.  This 
format  of  data  translation  is  implemented  to  illustrate  in  each  image,  the  instantaneous 
values  of  the  parameters  deemed  important  in  the  breakup  process. 
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Figure  8.  Timing  Sequence  for  High  Speed  Cinematography 
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4.7  Test  Matrix 


Experimental  modeling  of  the  atomization  process  inside  a  high  pressure  liquid 
rocket  combustion  chamber  is  difficult  due  to  the  extreme  temperature  and  turbulent 
conditions  which  exist.  Under  noncombustion  conditions  the  best  an  experiment  can  do  is 
match  certain  fluid  properties  and  dimensionless  parameters  deemed  important  in  the 
atomization  process.  For  the  simulant  liquid  study,  the  following  parameters  are  chosen 
based  on  their  importance  and  their  ability  to  match  hot-fire  conditions: 


1. 

2. 

3. 

4. 

5. 

6. 

7. 

8. 


Liquid  surface  tension  (a) 

Liquid  viscosity  (vj) 

Liquid  to  chamber-gas  density  ratio  (pi/pcg) 
Aerodynamic  Weber  number  (pgDVg^/a) 

Gaseous  Reynolds  number  (pgDVg/jxg) 

Liquid  to  gas  dynamic  pressure  ratio  (piVj^/pgVg^) 
Liquid  to  gas  kinematic  viscosity  ratio  (V\/Vg) 
Liquid  jet  diameter  (D) 


A  test  matrix  was  developed  to  identify  shock  tube/injector  operating  conditions 
which  simultaneously  match  most,  if  not  all,  of  these  parameters  in  the  liquid  rocket  engine 
under  study. 


The  liquid  rocket  engine  (LRE)  being  simulated  is  the  liquid  oxygen/hydrogen 
(LOX/H2)  Space  Shuttle  Main  Engine  (SSME)  with  injectable  LOX  and  the  combustion 
products  as  the  primary  fluids.  Table  3  lists  fluid  properties  and  dimensionless  parameters 
present  in  this  engine.  The  expected  chamber  properties  have  been  defined  so  that  proper 
cold  flow  studies  are  conceived.  The  cold  flow  properties  are  set  up  such  that  the  injected 
conditions  and  the  chamber  conditions  are  chosen  to  match  the  parameters  of  the 
combustion  products  in  the  SSME  as  well  as  possible.  To  effectively  do  so,  a  certain 
range  of  shock  strengths  have  initially  been  assumed.  With  this  information  all  relevant 
gas  properties  have  been  obtained  by  using  standard  shock  tube  relations  (Appendix  A). 
For  each  liquid  simulant  its  associated  liquid  properties  (e.g.  surface  tension  and 
viscosity)  have  been  used  to  determine  Weg,  Reg,  liquid/gas  viscosity,  density,  and 
dynamic  pressure  ratio  for  each  shock  tube  condition.  The  parameter  values  obtained  are 
then  compared  to  the  ones  present  in  the  engine  under  study.  The  parameter  values  in 
range  of  those  found  in  the  LOX/H2  engine  are  then  selected  along  with  its  associated 
initial  gas  pressure  of  the  driver  and  the  driven  section.  Thus,  denoting  operating 
conditions  which  will  match  the  aforementioned  parameters  in  this  LOX/H2  engine. 
Details  of  the  test  matrix  calculations  are  included  in  Appendix  A.  Various  wave 
strength  and  duration,  jet  velocity  and  diameter,  and  LOX  simulants  have  been  used  in  the 
delineation  of  the  breakup  mechanism  at  subcritical  conditions. 

A  cylindrical  driver  section  was  used  to  produce  a  step  wave  of  long  gas  flow 
duration,  such  that  the  pressure  behind  the  wave  was  relatively  constant  until  the 
reflected  shock  obstructed  the  flow.  Thus,  the  liquid  jet  underwent  positive  acceleration 
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during  the  time  that  the  images  were  taken.  In  addition,  images  of  jet  breakup 
undergoing  both  positive  and  negative  accelerations  were  also  taken.  This  type  of 
motion  is  accomplished  by  exposing  the  liquid  jet  to  N-waves  (Temkin,  1982).  In  these, 
the  fluid  velocity  decreases  from  some  value  immediately  behind  the  leading  shock,  to  a 
value  of  similar  magnitude,  but  having  the  opposite  direction,  immediately  ahead  of  the 
trailing  the  shock.  Thus,  owing  to  the  passage  of  the  wave,  the  liquid  jet  will  first 
experience  large  positive  accelerations  which  impart  to  the  jet  a  velocity  in  the  direction 
of  the  flow,  producing  initially  a  decrease  of  the  relative  fluid  velocity.  However,  owing 
to  the  deceleration  of  the  fluid  velocity  in  the  wave,  the  relative  velocity  will  eventually 
increase  in  time. 


5.0  RESULTS 


5.1  Overview 

Previous  studies  of  jet  breakup  in  a  crossflow  have  been  limited  to  low  time 
resolution  images  on  the  order  of  milliseconds.  Missing  are  the  details  of  the  breakup 
process,  in  general,  and  the  incipient  processes  in  particular.  In  addition,  studies  have  not 
been  conducted  that  study  the  breakup  of  a  liquid  jet  in  an  exponentially  decaying  velocity 
crossflow  (N-wave).  In  past  studies,  breakup  has  been  addressed  only  in  a  constant  gas 
velocity  flow  (step  wave).  Finally,  in  past  experiments  where  photographic  images  of  the 
atomization  process  have  been  taken,  only  large  liquid  to  gas  momentum  ratios  (q  >  1) 
have  been  addressed.  In  this  study,  the  simulation  of  a  liquid  rocket  engine  operating 
condition  was  of  paramount  concern.  In  order  to  do  so,  low  value  of  momentum  ratios 
were  obtained  (q  «  1).  This  section  presents  the  results  of  a  qualitative  and  quantitative 
study  that  addressed  the  effect  of  shock  wave  type  and  strength,  and  jet  diameter  on  jet 
breakup.  Shadowgraph  images  of  33  microsecond  (|is)  time  resolution  and  20 
nanosecond  (ns)  exposure  time  were  acquired  and  used  to  obtain  the  necessary 
information. 

A  test  matrix  was  developed  and  used  to  establish  operating  conditions  that  best 
match  the  parameters  deemed  important  in  the  breakup  process  present  in  the  SSME.  In  a 
liquid  rocket  engine,  very  high  Weber  and  Reynolds  numbers  are  present  due  to  the  highly 
turbulent  flow  conditions  and  the  fluid  properties  of  the  liquid  oxidizer  (low  surface 
tension).  In  the  experiments  conducted  in  the  laboratory,  very  high  speed  gas  flow 
velocities  were  difficult  to  obtain  due  to  hardware  limitations.  In  addition,  low  surface 
tension  simulants  were  found  to  be  challenging  in  their  storage  and  their  injection  into  the 
shock  tube  test  section.  As  a  result,  the  operating  conditions  selected  encompassed  a 
range  of  Weber  and  Reynolds  numbers  that  is  lower  than  the  SSME  value,  but  within  one 
order  of  magnitude. 
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Table  -  3  SSME  L0X/H2  Injectable  LOX  -  Combustion  Products 


LIQUID  OXYGEN  (LOX)  PROPERTIES 


Surface  tension 


Dynamic  viscosity  (ke/ms' 


Injection  Temperature  (K' 


GAS  PROPERTIES 


Gas 


Shock  strength 


Gas  flow  velocity  (m/s 


Chamber  density  (kg/m3 


Dynamic  viscosity  (kg/m  s 


Specific  heat  ratio 


Molecular  weight 


Speed  of  sound  (m/s 


Density  behind  shock  (kg/m3 


Shock  Mach  number 


FLUID  PROPERTIES 


Liquid  Reynolds  number 


Gaseous  Reynolds  number  (pgDVg/Up) 


Aerodynamic  Weber  number  (pgDVp^/a 


Liquid/chamber  gas  density  ratio  (pi/prg) 


Ohnesorge  number 


Liquid/gas  yiscosity  ratio 


Liquid/gas  dynamic  pressure  ratio 


Liquid/gas  kinematic  viscosity  ratio  (vi/v 


ulmSnB 


Injection  velocity  (m/s' 

) 

29.47 

Injector  diameter  (mm) 

1 

4.78 

1153 


0.0121 


0.00022 


111 


Combustion 

Products 


1.1  -2.0 


133  -  1000 


14.7 


1.1481 


1940 


13.589 


1595.6 


16  -  26.75 


1.05  - 1.40 


3.4375 


0.07  -  6.64 


0.0476  -  0.0797 


Figure  9  is  a  summary  of  the  maximum  Weber  -  Reynolds  number  plot  for  all  of 
the  conditions  tested.  The  Weber  -  Reynolds  number  time  history  for  each  individual  test 
is  illustrated  in  Appendix  C.  Figure  10  illustrates  the  minimum  dynamic  pressure  ratio  and 
maximum  viscosity  ratio  distribution  for  the  experiments.  The  shaded  box  denotes  the 
range  of  those  values  present  in  the  SSME.  The  hquid/gas  dynamic  pressure  ratio 
obtained  from  these  experiments  ranged  from  the  SSME  value  to  one  order  of  magnitude 
less.  However,  in  all  of  the  experiments  conducted,  the  liquid/gas  kinematic  viscosity  ratio 
was  always  within  the  range  present  in  the  SSME.  A  summary  of  the  test  conditions  is 
presented  in  the  Table  4. 

Instantaneous  pressure  measurements  were  collected  and  used  to  determine 
properties  of  the  gas  flow,  such  as,  velocity,  density,  and  temperature  dependent- viscosity. 
Values  for  the  aforementioned  parameters,  for  example,  were  calculated  using  the 
pressure  plot  from  the  pressure  transducer  upstream  of  the  hquid  injector.  These  values 
were  obtained  from  normal-shock  relations.  The  dynamic  viscosity  of  the  gas  was 
calculated  from  temperature  variations  across  the  shock.  From  this  the  following 
dimensionless  parameters  were  calculated; 


Weber  Number; 


Reynolds  Number; 


Liquid/Gas  Dynamic  Pressure  Ratio; 


Liquid/Gas  Kinematic  Viscosity  Ratio; 


We  = 


Re  = 


a 

Q  V  D. 

8  J 


^  p 


VR  = 


P/  P. 

P*  Pi 


(5-1) 

(5-2) 

(5-3) 

(5-4) 


The  Weber  number  characterizes  the  surface-tension  effect.  It  is  the  ratio  of  the 
dynamic  force  of  the  gas  to  the  surface-tension  of  the  hquid.  The  Reynolds  number  is  the 
ratio  of  the  inertial  force  to  the  viscous  forces  of  the  gas.  An  important  parameter  related 
to  the  behavior  of  the  jet  column  is  the  dynamic  pressure  ratio  which  is  simply  the 
momentum  flux  ratios  between  the  two  fluids.  Finally,  the  viscosity  ratio,  describes  the 
viscous  diffusivity  between  the  two  fluids,  and  is  composed  entirely  of  fluid  properties,  not 
geometric  or  flow  parameters.  The  ratios  are  dimensionless  fluid  properties  which  give  a 
measure  of  relative  rate  of  diffusion. 
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Teste  ETest7  ®Test8  ©  Test  9 


Pressure  Ratio  (q)  vs  Viscosity  Ratio  Distribution 
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The  behavior  of  the  jet  was  examined  under  low  liquid  to  gas  momentum  flows  in 
both  a  constant  and  an  exponentially  decaying  flow  field.  Highly  time  resolved  images  of 
the  primary  breakup  of  a  liquid  jet  were  obtained  to  provide  information  regarding  the 
different  regimes  a  jet  undergoes  before  breakup.  In  addition,  still  images  were  obtained 
in  Tests  1  and  2  and  used  to  obtain  drop  size  measurements.  The  rotating  drum  camera 
was  used  for  Tests  3  -  9.  The  time  duration  between  each  frame  was  33  |is  with  an 
exposure  time  of  20  nanoseconds.  Qualitative  and  quantitative  data  analyses  were 
employed  to  investigate  the  following  characteristics  of  the  breakup: 

•  Distortion  of  liquid  cylinder 

•  Jet  displacement 

•  Ligament  and  drop  formation 

•  Stripping  of  small  drops 

•  Surface  wave  formation 

•  Spreading  angle 

•  Jet  breakup  time 

•  Droplet  size  and  mass  distribution 

The  results  are  presented  in  the  following  sequence: 

•  Shock  wave  type  (Square  or  N-wave) 

•  Shock  strength  (from  34%  to  254%) 

•  Jet  diameter  (0.8  mm  to  1.41  mm) 

•  Surface  wave  formation 


5.2  Square  Wave  Perturbation 

The  square  wave  produced  a  relatively  constant  pressure.  As  a  result,  a  relatively 
constant  gas  velocity  was  realized  during  the  time  frame  in  which  the  images  were 
captured.  Figures  11-14  illustrate  the  jet  breakup  sequence  for  square  wave  Tests  6-9, 
respectively.  The  results  reveal  that  ligament  formation  occurs  only  after  that  jet  has 
distorted  to  an  elliptical  shape.  As  the  jet  continues  to  deform,  a  liquid  sheet  is  stripped 
away  from  the  outer  edge  of  the  liquid  core.  Only  after  the  jet  is  completely  flattened  to  a 
liquid  sheet  does  displacement  from  the  injector  occur.  At  high  Weber  numbers,  this 
observation  can  be  modified.  Here,  surface  waves  are  formed  along  the  sheet  plane 
which  induce  more  breakup  further  downstream.  During  the  entire  ligament/sheet 
formation  process,  drops  are  stripped  away  and  are  entrained  in  the  gas  flow.  Average 
size  of  drops  are  30  -  50  |jm. 


5.2.1  Jet  Displacement 

A  plot  of  the  data  obtained  from  experiments  with  a  square  wave  in  the 


nondimensional  form  of  X  =  x/Di  as  a  function  of  7"  —  ——  . 

Experiments  with  square  waves  have  a  constant  gas  flow  velocity  immediately  behind  the 


I—  is  shown  in  Figure  15. 

P. 
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shock  front.  The  displacement  is  a  result  of  a  normal  stress  acting  on  the  jet  column 
caused  by  the  pressure  drop  immediately  behind  the  wave.  The  values  of  the  downstream 
displacement  of  the  jet  element  as  a  function  of  its  time  after  injection  y/Uj,  where  x  and  y 
are  the  coordinates  of  the  jet  measured  from  the  orifice  parallel  and  normal  to  the 
freestream  velocity,  respectively,  are  discussed  and  compared  with  N-wave  experiments  in 
Section  5.3.1. 


5.2.2  Jet  Breakup  Time 

The  dimensionless  time  for  jet  breakup  was  derived  in  Section  2.5  (Reinecke, 
1978).  The  value  of  K  was  evaluated  for  experiments  with  Weber  number  over  600  and 
constant  gas  flow  velocity.  The  value  of  K  is  derived  from  knowledge  of  gas  flow 
velocity  and  density,  time  for  jet  shattering,  and  liquid  jet  diameter  and  density.  The  value 
of  K  could  not  be  evaluated  for  experiments  with  a  decaying  gas  flow  velocity  (N-waves), 
because  K  is  no  longer  a  constant.  For  the  square-wave  experiments  the  following 
observations  were  made: 


Test 

K 

Tj 

Weber  Number 

6 

9.22 

7.72 

907 

7 

5.023 

4.20 

3290 

8 

4.50 

3.76 

701 

9 

7.035 

5.89 

10400 

The  average  value  of  K  and  Tj  were  6.44  and  5.40,  respectively.  Results  from 
experiments  conducted  by  Reinecke  (1978)  were  similar  (K  =  7.0  and  Tj  =  5.85). 

5.3  N-Wave  Perturbation 

A  series  of  experiments  were  conducted  using  a  conical  driver  section  to  produce  a 
steep-fronted  N-wave  which  was  found  to  best  represent  an  actual  instability  wave  in  a 
liquid  rocket  engine.  A  36”  long  driver  section  was  used  to  produce  a  5.5  ms  wave 
duration.  A  10”  driver  was  used  to  generate  a  2.3  msec  wave  duration.  For  these 
experiments,  a  smaller  mid-diaphragm  holder  had  to  be  fabricated  to  enable  the  shock 
initiation.  The  conical  driver  section  volume  was  considerably  smaller  compared  with  the 
cylindrical  driver  section.  This  reduced  volume  had  no  effect  in  rupturing  the  second 
diaphragm  when  the  driver  gas  was  exhausted  into  the  intermediate  section  due  to  the 
reduction  of  pressure.  Thus,  a  smaller  mid-diaphragm  holder  increased  the  stored 
pressure  and  built  a  higher  pressure  differential  across  the  second  diaphragm  causing  it  to 
burst. 
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Figure  13.  Jet  breakup  sequence  for  Test  8 
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Figure  15.  Dimensionless  downstream  displacement  of  a  liquid  jet  under  the  influence  of  a  square  wave. 


Figures  16-20  illustrate  the  jet  breakup  sequence  for  N-wave  Tests  1-5, 
respectively.  The  results  from  these  experiments  indicate  that  the  atomization  controlling 
step  is  identical  to  that  behind  the  square  wave.  The  same  physics,  or  breakup  regime, 
apphes  to  both  square  and  N-wave  induced  breakup.  This  suggests  that  the  breakup  is 
controlled  by  the  initial  acceleration  of  the  gas  flow,  not  the  rate  of  decay. 


5.3,1  Jet  Displacement 


The  values  of  the  downstream  displacement  of  the  jet  element  as  a  function  of  its 
time  after  injection  y/Uj,  where  x  and  y  are  the  coordinates  of  the  jet  measured  from  the 
orifice  parallel  and  normal  to  the  freestream  velocity,  respectively,  are  available  and 
discussed  below  for  each  experiment.  A  plot  of  the  data  obtained  from  experiments  with 


an  N-wave  in  the  nondimensional  form  of  X  =  x/Dj  as  a  function  of 


IS 


shown  in  Figure  21.  Experiments  with  N- waves  have  a  decaying  gas  flow  velocity 
immediately  behind  the  shock  front.  The  liquid  jet  was  shown  to  displace  at  a  slower 
rate,  as  compared  with  its  displacement  behind  a  square  wave.  This  is  a  result  of  a 
decaying  normal  stress  acting  on  the  jet  column  caused  by  the  pressure  drop  immediately 
behind  the  wave. 


In  Tests  3  and  8,  the  same  shock  wave  strength  was  generated  (38%)  for  the  same 
jet  diameter  (1.04  mm)  but  with  a  2.3  ms  duration  N-wave  and  a  long  duration  square 
wave,  respectively.  The  percentage  denotes  the  percent  of  pressure  increase  behind  the 
shock  front.  In  other  words: 


Shock  Strength  (%)  = 


^AP 


xlOO 

J 


or 


Shock  Strength  (%)  = 


XlOO 


where  P2/P1  is  the  shock  pressure  ratio,  and  AP  is  the  pressure  difference  across  the 
shock. 


Pressure  plots  for  these  two  waves  are  shown  in  Figures  22  and  23,  respectively. 
The  observation  time  for  Test  3  and  Test  8  was  1.1  ms.  According  to  the  jet  breakup 
sequence  of  images  in  Figures  13  and  18  and  the  jet  displacement  plot  in  Figure  24,  there 
is  evidence  that  the  jet  responds  differently  to  the  N-wave.  This  is  supported  by  the  fact 
that,  in  the  displacement  plot,  the  jet  displacement  rate  is  less  for  the  N-wave  than  it  is  for 
the  constant  flow  induced  by  the  square  wave.  A  curve  fit  test  was  completed  on  the 
square  wave  plot  and  found  that  X  varied  with  T^  ^^  (Appendix  B).  With  a  constant  gas 
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Initial  Jet  t  -  60  ^is 


t  =  30  |is 


t  =  1 80  [IS 


Figure  16a.  Jet  breakup  sequence  for  Test  1 


35 


t  =  240  )j.s 


t  =  420  |is 


t  =  270  |is 


t  =  480  (IS 


t  =  360  (IS  t  =  540  (IS 

Figure  16B.  Still  images  of  different  stages  of  jet  breakup  for  Test  1 
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Figure  17a.  Jet  breakup  sequence  for  Test  2 
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Figure  17B.  Still  images  of  different  stages  of  jet  breakup  for  Test  2 
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Figure  18.  Jet  breakup  sequence  for  Test  3 
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Figure  19.  Jet  breakup  sequence  for  Test  4 
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Figure  20.  Jet  breakup  sequence  for  Test  5 
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Dimensionless  downstream  displacement  of  a  liquid  jet  under  the  influence  of  an  N-wave. 


flow  velocity  and  density,  the  value  of  T  will  only  increase  with  time  at  a  constant  rate. 
However,  the  value  of  X  is  changing  more  rapidly  denoting  that  there  is  an  acceleration  of 
the  jet  into  the  freestream  flow.  A  curve  fit  for  Test  3  resulted  in  X  varying  with  Ti-’o. 
Therefore  the  jet  accelerates  with  the  gas  flow  behind  an  N-wave  relative  to  the 
acceleration  experienced  with  a  square  wave. 


In  Figure  25,  the  dimensionless  displacement  is  plotted  for  another  condition 
where  the  N-wave  can  be  compared  to  the  square  wave  (Tests  4  and  6).  The  pressure 
plots  are  shown  in  Figures  26  and  27.  These  two  tests  were  similar  to  the  one  studied 
previously,  however,  a  larger  jet  diameter  was  used  (1.41  mm).  The  larger  jet  diameter 
appears  to  have  slowed  the  jet  displacement.  It  is  more  evident  here  that  the  N-wave 
decreased  the  value  of  T  due  to  the  sudden  drop  of  gas  flow  velocity  and  density.  In  the 


expression  for  T  =  ,  the  value  of  the  product  of  the  gas 

root  of  the  gas  density,  decreases  at  a  faster  rate  behind 


velocity  and  the  square 
the  decaying  flowfield 


on  an  N-wave  compared  to  the  rate  of  increase  of  t.  Meanwhile,  the  jet  continues  to  get 
displaced.  As  a  result,  dX/dT,  in  Figure  25,  will  initially  have  a  positive  slope  and  then 
have  a  negative  slope  when  the  gas  velocity  behind  the  wave  decays  rapidly.  A  similar 
study  was  conducted  between  Tests  5  and  7  but  at  a  higher  shock  strength  (Figures  28  and 
29,  respectively).  The  profile  of  the  displacement  for  the  square  wave  is  similar  to  that 
observed  at  the  lower  gas  velocity  (Figure  30).  However,  the  N-wave  profile  differs.  The 
gas  flow  behind  the  N-wave  is  higher  than  it  was  for  Test  4.  Therefore,  the  jet  displaced 
quicker  and  limited  the  observation  time.  During  this  time  the  N-wave  did  not  have  a  very 
strong  effect  in  decelerating  the  jet  displacement  Therefore,  the  gas  flow  velocity 
remained  almost  constant.  This  explains  the  reason  the  two  profiles  are  rather  similar  to 
each  other.  In  other  words,  during  the  observation  time.  Test  5  behaved  similar  to  Test  7 
since,  during  the  limited  observation  time,  the  gas  velocity  did  not  significantly  decline 
behind  the  N-wave. 


5.3.2  Ligament  and  Drop  Formation 

A  detailed  analysis  of  ligament  and  drop  formation  was  undertaken  for  Tests  1  and 
2.  These  images  were  captured  using  high  resolution  still  photographs,  and  were  visually 
inspected  for  the  different  regimes  of  breakup  (Figures  16  and  17).  Table  4  illustrates  the 
gas  and  liquid  phase  properties  calculated  for  these  tests.  Test  1,  with  a  peak  pressure 
increase  of  34%,  had  a  shock  Mach  number  of  1.14  with  a  peak  gas  flow  velocity  of  76 
m/s.  The  pressure  time  history  is  shown  in  Figure  31.  The  pressure  plots  also  represent 
the  time  and  the  instantaneous  pressure  value  for  each  picture. 
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Figure  22.  Pressure  Time  History  for  Test  3  (top)  full  duration  (bottom)  1.5  msec  duration 
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Figure  23.  Pressure  Time  History  for  Test  8  (top)  full  duration  (bottom)  1.5  msec  duration 
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8  S  38  infinite  1.04 


Figure  24.  Dimensionless  downstream  displacement  of  a  liquid  jet  for  Test  3  and  Test  8. 
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Figure  25.  Dimensionless  downstream  displacement  of  a  liquid  jet  for  Test  4  and  Test  6 
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Figure  27.  Pressure  Time  History  for  Test  6  (top)  full  duration  (bottom)  1.5  msec  duration 
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The  methanol  jet  diameter  was  1.0  mm  with  a  jet  velocity  of  1.0  m/s.  Figures  16a 
and  16b  illustrate  stUl  images  of  the  breakup  sequence  for  Test  1.  Note,  that  each  image 
for  Test  1  and  2  is  a  result  of  a  different  experiment  captured  with  a  35  mm  Pentax  camera 
at  various  stages  of  breakup.  Each  experiment  produced  repeatable  shocks  to  within  0.1 
psi. 


In  this  subsonic  crossflow,  the  jet  expands  uniformly  and  breaks  off  at  the  trough 
of  a  column  wave.  The  images  show  that,  with  the  passage  of  an  N-wave,  the  jet 
undergoes  distortion  during  the  0  -  60  |is  time  frame.  Ligament  formation  occurs  at  the 
outer  edge  of  the  ellipsoid  liquid  core  at  120  \is.  Thin  liquid  sheets  become  visible  at  this 
time  at  which  point  50  fim  drops  are  stripped  away  and  entrained  in  the  gas  flow.  After 
360  |is,  surface  waves  appear  on  the  liquid  sheet  surface  due  to  the  turbulent  gas  flow,  and 
jet  displacement  is  initiated.  Clouds  of  drops  are  entrained  in  the  gas  flow  as  the  surface 
wave  amplitude  increases.  The  remainder  of  the  jet  is  entrained  in  the  gas  flow  with  50  - 
75  ^im  droplets.  The  last  image,  taken  at  540  \is,  did  not  encompass  the  complete  breakup 
process  and  the  behavior  of  the  jet  in  a  decaying  gas  flow  velocity  field.  In  fact,  the  gas 
velocity  remained  relatively  constant  during  these  times. 

Test  2  was  conducted  with  a  higher  initial  diaphragm  pressure  ratio  and  resulted  in 
a  57%  peak  pressure  increase  across  the  shock  wave.  The  shock  Mach  number  was  1.22 
with  peak  gas  flow  velocity  of  122  m/s.  Figure  32  shows  the  pressure  time  history  for  the 
57%  wave.  The  jet  breakup  sequence  is  shown  in  Figure  17a  and  17b.  At  30  p,s  after  the 
passage  of  the  shock  wave,  the  jet  initiates  the  distortion  process  and  flattens  into  an 
ellipse  shape.  Thin  ligaments  form  at  the  outer  edge  of  the  ellipsoid  after  60  [is.  Drops 
begin  to  entrain  in  the  gas  flow  at  90  its.  After  120  [is,  the  liquid  core  has  completely 
flattened  to  a  long  liquid  sheet.  The  average  size  of  drops  entrained  in  the  gas  flow  is  -30 
|im  The  jet  displaces  from  the  injector  after  180  p-s  and  travels  downstream  of  the  shock 
tube  at  16  m/s.  Surface  waves  form  on  the  sheet  surface  at  this  time.  A  cloud  of  fine 
droplets  are  stripped  away  from  the  liquid  sheet  and  entrained  in  the  gas  flow.  At  270  |is, 
the  liquid  sheet  has  displaced  3  mm  away  from  the  injector  centerline. 

After  540  |is,  fresh  liquid  has  emerged  from  the  injector  exit.  This  indicates  that 
the  gas  flow  behind  the  shock  has  no  effect  on  jet  recovery  rate,  which  was  approximately 
equal  to  the  nominal  jet  velocity.  The  shock  front  may,  however,  produce  an  interruption 
of  the  liquid  flow  into  the  test  section. 

Droplet  diameters  were  measured  for  a  selected  area  of  some  of  the  jet 
photographs  for  Tests  3, 4, 5, 6, 7,  8,  and  9.  The  droplets  were  measured  by  enlarging  the 
selected  area  which  results  in  a  total  magnification  of  70,  and  measuring  the  diameters 
with  a  graduated  scale.  The  droplet  diameter  determination  obtained  here  has  validity  only 
as  a  determination  of  the  sizes  of  initial  droplets  (i.e.,  droplets  close  to  the  point  of 
injection  as  is  the  case  for  all  of  the  particles  in  the  photographs)  for  each  test.  The  degree 
of  jet  breakup  was  not  a  consideration  in  these  measurements,  and  may  be  very  large  or 
very  small.  Furthermore,  the  downstream  droplet  distribution  would  probably  be  quite 
different,  due  to  secondary  breakup  and  evaporation.  The  measurement  was  limited  to 
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Figure  32.  Pressure  Time  History  for  Test  2  a)  Full  wave  duration 
b)  1  msec  wave  duration 


particles  which  approximate  spheres.  This  was  arbitrarily  defined  as  particles  for  which  the 
maximum  diameter  does  not  exceed  the  minimum  diameter  by  more  than  two  to  one.  For 
the  nonspherical  particles  thus  counted,  the  droplet  diameter  was  defined  by  (max-diam  + 
min-diam)/2.  Furthermore,  only  those  particles  and  droplets  that  were  in  good  focus  were 
counted. 

Under  the  above  restrictions,  the  mean  droplet  diameters  were  measured  and 
tabulated  in  Table  5.  Figure  33  illustrates  the  beneficial  influence  of  increased  levels  of 
airstream  velocity  upon  the  spray  drop  size  for  two  injector  diameters.  Figure  34 
illustrates  the  influence  that  the  diameter  of  the  injection  orifice  has  upon  the  spray  quality. 
The  average  dropsize  is  observed  to  increase  with  an  increase  in  orifice  diameter. 


Table  5.  Mean  measured  droplet  diameters 


Test  Number 

Jet  Diam.  (mm) 

Jet  Veloc.(m/s) 

Gas  VeIoc.(m/s) 

1 

30 

0.8 

1.0 

76 

2 

22 

1.0 

1.0 

122 

3 

37 

1.04 

0.71 

89 

4 

50 

1.41 

0.435 

85 

5 

33 

1.41 

0.435 

135 

6 

43 

1.41 

0.435 

97 

7 

28 

1.41 

0.435 

168 

8 

26 

1.04 

0.71 

99 

9 

17 

1.04 

0.71 

300 

5.4  Parametric  Studies 

In  addition  to  exploring  two  wave  types,  variations  in  the  following  parameters 
were  explored: 

•  Shock  wave  strength 

•  Jet  diameter 

5.4.1  Effect  of  Shock  Strength 


a  higher  shock  strength  produces  a  more  rapid  breakup  than  the  34%  wave.  In  addition, 
quantitative  analysis  was  done  by  plotting  the  jet  displacement  data  in  nondimensional 


form 


of  Z  =  x/Dj  as  a  function  of 


D; 


Since  the  images  were  all  captured 


before  the  jet  was  influenced  by  the  negative  acceleration,  (the  last  data  point  for  Tests  1 
and  2  were  taken  at  t=540  |is)  the  assumption  is  made  that  the  freestream  velocity  was 
relatively  constant.  Based  on  this  assumption,  the  displacement  plot  for  two  conditions 
was  expected  to  have  a  similar  slope.  However,  the  displacement  of  the  jet  behind  a  57% 
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Figure  33.  Influence  of  airstream  velocity  upon  spray  dropsize 
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Figure  34.  Influence  of  Liquid  Orifice  Diameter  Upon  Spray  Dropsize 


wave  occurred  much  quicker.  According  to  the  displacement  plot,  the  57%  wave 
produced  a  much  higher  displacement  value,  at  the  time  the  last  picture  was  taken.  This 
can  be  explained  by  the  dynamic  force  of  the  gas  flow  which  is  much  higher  for  the 
higher  strength  shock.  And  since  the  liquid  jet  dynamic  pressure  was  the  same  for  the 
two  cases,  the  momentum  transfer  between  the  gas  and  liquid  is  much  greater  for  the 
higher  strength  shock.  Therefore,  the  liquid/gas  dynamic  pressure  ratio  (q)  is  much 
smaller  for  the  57%  wave.  Hence,  displacement  occurs  much  more  rapidly.  The 
maximum  Weber  numbers  for  Tests  1  and  2  were  370  and  870,  respectively.  However, 
the  Weber  number  has  shown  to  have  greater  effect  on  jet  breakup  and  almost  none  on  jet 
displacement. 

The  effect  of  wave  strength  on  jet  breakup  was  also  studied  for  an  N-wave,  but  for 
a  different  wave  duration.  Tests  4  and  5  generated  a  2.3  miUisecond  (ms)  wave  duration 
with  a  37%  and  62%  shock  strength,  respectively.  The  pressure  plots  are  shown  in 
Figures  36  and  37.  Due  to  the  short  wave  duration  the  pressure  decrease  behind  this  type 
of  an  N-wave  is  much  more  sudden.  The  diameter  of  the  jet  was  1.41  mm.  From  the 
images  (Figures  19  and  20)  and  the  displacement  plot  in  Figure  38,  it  was  shown  that  the 
rate  of  jet  displacement  was  much  higher  for  the  62%  N-wave.  Due  to  this  faster 
displacement,  the  observation  time  of  the  jet  was  less  than  for  the  37%  wave  (725  p-s  vs 
1500  ps  ).  Therefore,  during  this  time  the  liquid  jet  did  not  have  the  opportunity  to 
decelerate  behind  the  N-wave.  In  fact,  there  was  only  a  20%  reduction  in  static  pressure 
of  the  gas  flow  at  725  ps  whereas,  for  the  37%  N-wave,  the  pressure  decreased  by  75%, 
denoting  a  large  deceleration  of  the  jet  by  1500  ps.  Due  to  this  pressure  drop,  the  gas 
velocity  and  density  also  experience  a  sudden  decrease. 

Since  the  dimensionless  time  T  is  dependent  on  this  decreasing  gas  velocity  and 
density,  the  value  of  T  begins  to  decrease  with  time  as  the  jet  displacement  continues  but 
at  a  slower  rate.  In  Figure  38,  this  behavior  is  clearly  illustrated  as  a  sudden  increase  of 
the  slope  of  the  plot  for  Test  4.  In  Test  5,  the  gas  flow  velocity  remained  relatively 
constant;  therefore,  with  an  increase  in  time  the  total  value  of  T  increased,  even  though 
the  gas  velocity  and  density  may  have  decreased  by  a  small  value.  Jet  behavior  was  further 
illustrated  in  the  comparison  between  Tests  4  and  7,  where  in  Test  7  a  65%  square  wave 
was  generated  (Figure  39).  The  gas  flow  velocity  behind  the  square  wave  was  constant  so 
that  the  value  of  T  increased  with  time.  This  explains  the  reason  why  the  slope  of  the 
dimensionless  downstream  displacement  in  Figure  40  is  smaller  for  the  constant  velocity 
Test  7.  The  jet  breakup  sequence  for  Test  7  is  shown  in  Figure  12. 

The  effect  of  liquid  jet  displacement  behind  a  low  and  high  strength  square  wave 
was  also  analyzed.  In  Test  8,  a  39%  square  wave  was  generated  (Figure  41).  In  Test  9  a 
254%  square  wave  was  generated  (Figure  42).  The  displacement  plot  for  the  1.04  mm 
diameter  jet  for  both  tests  is  shown  in  Figure  43.  Due  to  the  very  high  gas  flow  velocity 
behind  the  254%  wave  (300  m/s)  jet  displacement  occurred  very  rapidly.  Therefore  the 
observation  time  was  limited  to  460  ps  before  the  image  exited  the  view  field  of  the 
camera.  From  the  jet  breakup  sequence  for  Test  9  in  Figure  14,  it  is  clearly  visible  that  the 
jet  displaces  more  rapidly  as  compared  with  Test  8  (Figure  13). 
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Figure  35.  Dimensionless  downstream  displacement  of  a  liquid  jet  for  Test  1  and  Test  2 


In  summary,  higher  strength  shock  waves  displace  the  liquid  jet  faster  than  the 
lower  strength  waves.  This  was  found  to  be  independent  of  wave  type.  However,  higher 
strength  shock  displaced  the  jet  so  rapidly  that  the  jet  did  not  respond  to  the  decelerating 
gas  flow.  In  lower  strength  shocks,  the  jet  displaced  more  slowly,  allowing  it  to 
decelerate  in  the  decaying  gas  flow. 


5.4.2  Effect  of  Jet  Diameter 

Tests  3  and  4  were  compared  with  each  other  to  study  the  effect  of  jet 
displacement  with  a  change  in  jet  diameter.  Tests  3  and  4  were  both  a  38%  N-wave  of  2.3 
ms  wave  duration  (Figures  22  and  26).  However,  Test  3  was  conducted  with  a  1.04  mm 
and  Test  4  was  with  a  1.41  mm  jet  diameter.  A  comparison  of  the  jet  breakup  sequence 
for  Tests  3  and  4  in  Figures  18  and  19,  respectively,  illustrates  that  the  larger  jet  diameter 
causes  a  slower  displacement  rate,  and  therefore  a  longer  observation  time  (by  400  |is). 
During  this  additional  time  the  jet  decelerated  more.  Figure  44  shows  the  displacement 
plot  for  Tests  3  and  4.  The  additional  decrease  of  pressure  in  Test  4  contributed  to  a 
lower  value  of  T  due  to  the  lower  velocity  and  density  of  the  gas  flow.  The  values  of  X 
and  T  decrease  with  larger  diameter  jets.  Therefore,  the  smaller  diameter  jet  (1.04  mm)  in 
Test  3  would  result  in  a  larger  value  of  X  and  T.  This  explains  the  "stretching"  of  the  plot 
for  Test  3.  In  experiments  with  square  waves  the  similar  observation  was  found.  In  the 
displacement  plots  for  Tests  6  and  8  (Appendix  B)  it  is  shown  that  the  smaller  jet  diameter 
wiU  displace  from  the  injector  quicker  than  a  larger  diameter  jet.  However,  in  these 
experiments  there  is  no  decay  in  gas  flow  velocity. 

In  summary,  the  larger  the  jet  diameter,  the  more  resistant  it  is  to  displacement.  As  a 
result,  larger  diameter  jets,  in  a  decaying  gas  flow,  respond  more  to  this  negative 
acceleration  than  smaller  diameter  jets. 


5.5  Surface  Wave  Formation 

Before  the  passage  of  a  shock,  the  liquid  jet  is  a  steady,  circular,  coherent  column. 
After  passage  of  the  shock  and  the  onset  of  a  high  momentum  gas  flow,  axial  waves 
develop  along  the  surface  of  the  jet  column  and  propagate  along  the  plume  with  increasing 
ampUtude.  The  plume  curves  downstream,  depending  on  the  liquid/gas  dynamic  pressure 
ratio  (due  to  the  aerodynamic  drag)  and  eventually  fractures  at  the  trough  of  a  high- 
amplitude  wave.  Aerodynamic  forces  quickly  decompose  the  liquid  fragment.  The  spray 
then  becomes  atomized  within  a  short  distance. 
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Figure  38.  Dimensionless  downstream  displacement  of  a  liquid  jet  for  Test  4  and  Test  5 


T 


Test  Wave  Type  Strength  (%)  Duration  (ms)  Diameter  (mm) 

4  N  37  2.3  1.41 

7  S  65  infinite  1.41 


Figure  40.  Dimensionless  downstream  displacement  of  a  liquid  jet  for  Test  4  and  Test  7. 
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Pressure  (psi) 


Figure  41.  Pressure  Time  History  for  Test  8  (top)  full  duration  (bottom)  1.5  msec  duration 
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Picture  Number 


Figure  42  .  Pressure  Time  History  for  Test  9  (top)  full  duration  (bottom)  0.5  msec  duration 
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T 

Test  Wave  Type  Strength  (%)  Duration  (ms)  Diameter  (mm) 

8  S  39  infinite  1.04 

9  S  254  infinite  1.04 


Figure  43.  Dimensionless  downstream  displacement  of  a  liquid  jet  for  Test  8  and  Test  9. 
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The  wavelength  and  maximum  amplitude  of  the  surface  waves  observed  for  all  the 
tests  conducted  are  documented  in  Table  4.  A  plot  of  surface  wave  amplitude  as  a 
function  of  liquid/gas  dynamic  pressure  ratio  for  Tests  4, 5,  6,  and  9  is  shown  in  Figure  45. 
The  surface  wave  amplitudes  for  Tests  4,  5,  6,  and  9  were  measured  from  the  jet  breakup 
images  in  Figures  19,  20,  11,  and  14,  respectively.  The  plot  illustrates  that,  with  an 
increase  in  dynamic  pressure  ratio,  the  amplitude  of  the  surface  waves  reduces,  for  q  < 
0.002.  These  results  are  similar  to  those  found  by  Schetz  and  Padhye  (1977)  who 
conducted  studies  at  lower  gas  momentum  flows.  They  also  found  that  the  wavelength  of 
surface  waves  also  reduces  with  increasing  dynamic  pressure  ratio.  However,  no  solid 
evidence  for  this  was  found  in  this  present  investigation. 


•  Test  4  o  Test  6  ■  Test  9  □  Test  5 


Figure  45.  Jet  surface  wave  amplitude  vs  liquid/gas  dynamic  pressure  ratio. 


6.0  SUMMARY,  AND  CONCLUSIONS 


6.1  Simunary 

A  fast-response  imaging  and  diagnostic  technique  has  been  utilized  to  investigate 
the  primary  mechanism  of  breakup  of  a  liquid  jet  behind  a  pressure  pulse  as  it  pertains  to 
the  evolution  and  development  of  combustion  instability  in  liquid  rocket  engines.  The 
study  utilized  high  speed  cinematography  to  delineate  the  primary  atomization  of  a  liquid 
jet  under  subcritical  conditions.  A  test  matrix  was  followed  that  explores  operating 
conditions  which  best  match  SSME  operating  conditions.  Time-sequenced  still  images 
verified  good  shock  repeatability. 
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The  results  show  that  the  flow  behind  a  shock  wave  moving  transverse  to  the 
liquid  jet  causes  disintegration  to  occur  simultaneously  at  all  points  along  its  length.  The 
liquid  initially  responds  to  the  high  speed  gas  flow  due  to  the  external  pressure  distribution 
surrounding  the  cylindrical  jet.  The  flow  around  the  jet  acts  like  flow  around  a  solid 
cylinder.  The  jet  then  flattens  and,  at  the  edges  of  the  flattened  jet,  drops  and  ligaments 
are  tom  off  by  the  combined  action  of  the  tangential  stress  and  surface  tension. 

From  the  images,  the  following  measurements  were  made:  the  time  required  for  jet 
breakup,  jet  displacement  from  the  injector,  and  the  droplet  diameter. 

Two  wave  types  were  explored:  A  square  wave  induced  breakup,  and  an  N-wave 
induced  breakup.  In  addition,  variation  in  the  following  parameters  was  explored: 


•  Shock  wave  strength 

•  Jet  diameter 


6.1.1  Square  Wave  Induced  Breakup 

The  square  wave  produced  a  relatively  constant  pressure  (hence  a  constant  gas 
velocity)  during  the  time  frame  the  images  were  captured.  Images  indicate  that  ligament 
formation  occurs  only  after  the  jet  has  distorted  to  an  ellipsoidal  shape.  As  the  jet 
continues  to  deform,  droplets  are  stripped  away  from  the  outer  edge  of  the  liquid  core. 
Once  the  jet  has  flattened  to  a  liquid  sheet,  displacement  from  the  injector  occurs.  Surface 
waves  are  formed  along  the  sheet  which  induce  subsequent  breakup  further  downstream. 
During  the  entire  ligamenl/sheet  formation  process,  drops  are  stripped  away  and  are 
entrained  in  the  gas  flow.  The  average  sizes  of  drops  for  Tests  1  and  2  were  30  to  50  pm. 


The  jet  displacement  behind  the  square  wave  was  nondimensionalized  as  a  function 


of  T  I—.  Comparison  of  this  result  with  Reinecke’s  (1978)  dimensionless  jet 

Dj  ipj 

displacement  data,  for  low  gas  momentum  flows  (Test  8),  reveals  no  similarity  to  the 
parabolic  least-square  regression  (X=0.8T^)  obtained  from  their  results.  The  difference  is 
attributed  to  the  fact  that  Reinecke’s  data  used  a  general  value  of  7.0  for  K,  in  Eq.  2.3,  for 
Weber  number  greater  than  600.  However,  Test  8  had  a  value  of  4.5  for  K  and  Weber 
number  of  700.  Although,  each  test  in  the  present  study  resulted  in  a  different 
displacement  plot,  the  plots  for  Tests  6,  7,  and  9  are  generally  similar  to  the  parabolic 
least-square  regression,  indicating  better  agreement  at  higher  Weber  numbers. 


From  the  nondimensional  value  of  jet  breakup  up,  typical  values  of  K  and  T  for 
Weber  numbers  above  600  are  6.40  and  5.40,  respectively.  These  values  are  similar  to 
those  found  by  Reinecke  (K  =  7.0  and  Tj  =  5.85)  for  Weber  numbers  above  500. 


The  effect  of  jet  diameter  on  displacement  was  also  analyzed.  It  was  found  that 
smaller  diameter  jets  accelerate  with  the  gas  flow  much  quicker  than  the  larger  diameter 
jets.  Furthermore,  this  rapid  displacement  decreases  the  observation  time  for  smaller 
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diameter  jets.  Hence,  the  displacement  of  the  jet  could  be  measured  only  before  the  gas 
flow  velocity  had  an  opportunity  to  decay.  In  some  cases,  the  gas  velocity  for  such 
experiments  was  somewhat  similar  to  those  for  equal  strength  square  waves.  Larger  jet 
diameters  had  a  higher  resistance  to  the  gas  flow.  Therefore,  the  displacement  rate  was 
smaller  and,  by  the  time  the  last  image  was  captured,  the  jet  experienced  the  decelerating 
gas  flow. 

The  surface  wave  formation  on  the  liquid  jet  surface  was  analyzed  as  well.  It  was 
found  that  with  an  increase  ia  hquid/gas  dynamic  pressure  ratio  {q),  the  surface  wave 
amplitude  decreased.  This  result  is  similar  to  that  found  by  Schetz  and  Padhye  (1977). 
They  observed  that  with  an  increase  in  q,  the  amplitude  and  wavelength  of  surface  waves 
were  reduced.  There  was  no  evidence  of  wavelength  reduction  with  increasing  q  in  this 
investigation,  which  evaluated  higher  gas  momentum  flows. 

The  general  atomization  characteristics  of  the  jet  were  observed  qualitatively.  It 
was  apparent  that  high  Weber  and  Reynolds  number  gas  flows  produced  finer  droplets  in 
the  atomization  process.  A  good  example  of  this  is  illustrated  in  the  breakup  sequence  for 
the  254%  wave  (Figure  14).  The  value  of  the  Weber  and  Reynolds  number  for  this  test 
are  an  one  order  of  magnitude  lower  than  that  found  in  the  SSME.  This  indicates  that  in 
actual  SSME  rocket  conditions,  jet  breakup  occurs  faster  than  that  depicted  in  Figure  14. 

6.1.2  N-Wave  Induced  Breakup 

The  breakup  images  indicate  that  the  atomization  controlling  step  is  identical  to 
that  behind  the  step  wave.  The  same  physics  applies  to  both  step  and  N-wave  induced 
breakup.  For  the  34%  wave  it  was  shown  that  ligaments  strip  away  from  the  outer  edge 
of  the  deformed  liquid  and  drops  are  later  produced  and  entrained  in  the  gas  flow.  The 
bulk  liquid  experienced  displacement  from  the  injector  after  the  liquid  flattened  to  a  thin 
sheet.  The  57%  wave  illustrated  a  much  faster  atomization  rate  than  the  34%  wave.  In 
the  case  of  ligament  and  sheet  formation  and  also  jet  displacement,  the  57%  wave  required 
half  the  time  the  34%  wave  required  to  accomplish  the  processes.  This  was  also  shown  in 
the  drift  rate  between  these  two  cases. 

The  result  of  the  nondimensional  displacement  of  the  liquid  jet  reveals  that  the 
liquid  jet,  depending  on  the  observation  time,  decelerates  with  the  deca)dng  gas  velocity 
and  density.  The  displacement  plots  for  N-waves  generally  show  a  steeper  slope  than  the 
square  wave  results.  This  is  attributed  to  the  fact  that  the  value  of  T  is  declining  with 
time,  due  to  the  decreasing  gas  velocity  and  density.  Hence,  a  steeper  slope  is  generated. 
Qualitative  analysis  from  the  jet  breakup  sequence  of  images  for  the  N-wave  and  square 
wave  experiments  delineate  this  lower  rate  of  displacement  for  the  N-wave. 

The  analysis  on  the  effect  of  jet  diameter  on  displacement  for  N-wave  experiments 
shows  similarity  with  the  results  obtained  for  the  square  wave  experiments;  in  that  smaller 
diameter  jets  accelerate  with  the  gas  flow  much  quicker  than  the  larger  diameter  jets. 
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The  atomization  of  the  jet  was  similar  to  that  found  for  the  square  waves.  Finer 
droplets  were  generated  in  high  Weber  and  Reynolds  number  gas  flows.  In  actual  rocket 
engine  conditions,  the  Weber  and  Reynolds  number  are  much  higher  than  the  values  in 
this  investigation.  This  implies  that  jet  breakup  occurs  very  rapidly  -  much  like  that 
found  in  Figure  38.  In  these  extreme  conditions,  the  jet  may  shatter  before  the 
deceleration  of  the  gas  flow  is  realized  behind  an  N-wave.  The  issue  then  becomes  “how 
do  the  drops  that  are  formed  from  the  jet  core  react  under  a  decelerating  flow  field?”. 

6.1.3  Drop  Size  Measurements 

Drop  diameters  were  measured  for  Test  1  and  2  by  analyzing  the  still  images. 
For  Tests  3-9  the  images  were  enlarged  by  a  factor  of  70  and  drop  sizes  were  measured 
with  a  scale.  The  droplet  diameter  determination  obtained  here  has  validity  only  as  a 
determination  of  the  sizes  of  initial  droplets  for  each  test.  For  two  different  injector 
diameters  it  was  found  that  the  average  drop  size  decreased  with  airstream  velocity.  The 
larger  diameter  jet  (1.41  mm)  produced  the  larger  drop  sizes  (on  the  order  of  50  fim). 
The  smallest  detected  drops  (15  pm)  were  found  for  the  case  where  the  gas  velocity  was 
300  m/s  and  the  jet  diameter  was  1.04  mm. 

6.2  Conclusions 

•  N-waves  decelerate  jet  displacement  and  the  jet  breakup  time  -  unless  very  rapid 
breakup  occurs. 

•  Faster  atomization  rates  and  jet  displacements  occurred  behind  a  step  wave. 

•  Jet  displacement  rate  increases  with  shock  wave  strength. 

•  Smaller  diameter  jets  displace  faster  from  the  injector  than  larger  diameter  jets. 

•  Surface  wave  amplitude  decreases  with  increasing  liquid/gas  dynamic  pressure  ratio. 

•  Drop  size  decreases  with  increasing  shock  wave  strength. 

•  Drop  size  increases  with  liquid  jet  diameter. 

•  Shock  interaction  with  the  primary  atomization  process  produces  a  substantial 
change  to  the  breakup  mechanism  and  serves,  as  a  result,  a  principal  candidate  for 
the  promotion  and  acceleration  of  rocket  engine  instability 


6.3  Reconunendations  for  Future  Work 

Although  these  experiments  provide  insight  on  the  breakup  of  a  liquid  jet,  several 
steps  can  be  taken  to  improve  the  understanding.  The  following  is  a  list  of 
recommendations  for  future  work; 


•  Conduct  experiments  at  elevated  pressures  in  the  driven  section  of  the  shock  tube. 

•  Liquids  with  lower  surface  tension  and  viscosity  should  be  used  to  more  accurately 
simulate  liquid  oxygen.  This  will  also  increase  the  Weber  number  and  Reynolds 
number  to  levels  closer  to  SSME  values. 
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•  Supercritical  simulant  liquids  should  be  used  to  evaluate  the  breakup  process  of  a 
liquid  jet  in  a  supercritical  environment.  This  can  be  done  by  injecting  the  liquid  at 
supercritical  pressures  but  subcritical  temperatures  into  the  shock  tube  which  is  at 
supercritical  pressure  and  temperature.  The  jet  will  then  experience  supercritical 
heating  as  it  penetrates  into  the  bore  of  the  shock  tube. 

•  Experiments  covering  the  entire  test  matrix  would  provide  the  most  useful 
information.  Specifically,  experiments  that  cover  the  entire  range  of  the  important 
fluid  parameters. 

•  A  laser  sheet  in  the  plane  of  the  jet  axis  can  capture  images  of  the  breakup  at  the 
focus  plane  of  the  camera.  Thus,  drops  that  are  formed  can  easily  be  seen  on  the 
image  plane.  Shadowgraphs  will  capture  images  of  drops  that  are  and  are  not  in 
focus,  making  drop  size  analysis  less  difficult. 

•  Scanning  of  the  image  and  the  implementation  of  an  image  analysis  software  will 
improve  the  analysis  of  the  breakup  and  drop  formation  -  provided  that  the  resolution 
of  the  scan  can  resolve  the  smallest  droplets. 
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Appendix  A:  Test  Matrix  Calculations 

1 .  Select  a  LOX  simulant  with  similar  surface  tension  and  viscosity. 

2.  Select  a  range  of  shock  strengths  (P2/P1). 

3.  For  each  shock  strength  value  calculate  gas  flow  velocity  behind  the  shock,  the 
shock  velocity,  gas  density  ratio  across  the  shock,  and  the  diaphragm  pressure  ratio 
using  the  equation  of  motion  for  a  normal  shock  moving  into  a  gas  at  rest  (John, 
1984) 


(A)  Moving 


V 

^2 

S 

•^2 

Pi 

(B)  Fixed 
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s  §2 
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^2 

The  equation  of  motion  for  a  normal  shock  moving  rightward  into  a  gas  at  rest 
(Figure  A)  is  written.  To  reduce  the  problem  to  one  of  steady  flow,  the  shock  is  fixed 
(Figure  B). 

It  follows  from  standard  shock  relations  that  the  gas  velocity  behind  the  shock  is 
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The  shock  velocity,  Vg,  is  solved  for  using  the  expression  in  A.1:. 


A.1 


V  =-iL(y  +  J)±  Li^ - 

^  4  ^  V  16 

It  follows  that  the  density  ratio  across  the  shock  is: 
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V  -V 
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A.3 


The  diaphragm  pressure  ratio  is  calculated,  considering  the  speed  of  sound  in  the 
driver  section  and  the  driven  section  of  the  shock  tube  is  the  same.  Setting  Vg2  = 
Vg3  and  P2  =  P3.  the  result  is: 
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4.  By  equating  the  We  equation  with  the  range  of  We  values  present  in  the  LRE  we 
determine,  for  each  shock  strength,  a  range  of  initial  gas  density  in  the  driven 
section.  Since  there  is  a  low  to  high  range  of  We  and  Vg2  and  q  are  fixed  values 
for  a  single  shock  strength  value  and  simulant,  respectively,  a  low  to  high  range  of 
P2  is  calculated.  This  is  denoted  on  the  spreadsheet  as  'Dl(low)’  and  'Dl(hi)’. 
Using  the  gas  density  ratio  calculated  in  step  3  the  gas  density  in  the  driven  section 
of  the  shock  tube,  pp  is  calculated. 

2 

Weber  =  =  5.96  X  ^  5.64  x  10^ 


Pi 


-9-(5.96x10'‘  ->5.64xl0^)*-r^ 

viD^  ^  fpil 


5.  Using  the  ideal  gas  law  equation  the  pressure  required  to  obtain  pi  in  the  driven 
section  is  then  obtained.  Again,  due  to  the  low  to  high  range  of  gas  densities  a 
range  of  low  to  high  driven  section  pressures  are  also  obtained. 


6. 


In  step  3  the  diaphragm  pressure  ratio,  (p4/pi),  was  calculated  for  each  shock 
strength  value.  Now  that  pj  is  known  from  the  previous  step,  P4  is  cdculated  from 


the  ratio  value.  A  low  to  high  range  of  p4  is  also  obtained. 


P4=Pl 


I 


1.  The  gaseous  Reynolds  Number  (Re)  can  be  calculated  using  the  gas  properties 
obtained  in  the  spreadsheet.  In  these  experiments  air  is  used  as  the  shock  tube  gas. 
A  low  to  high  range  of  Re  is  obtained  due  to  its  dependence  on  pg2. 

=  =  116x10^  -^1.46x10^ 

P2 

Behind  the  shock  wave  there  is  a  temperature  rise  in  the  gas.  This  temperature 
variation  causes  a  change  in  the  gas  viscosity.  Therefore,  the  viscosity  term  is  a 
function  of  shock  strength. 
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8.  The  liquid/chamber-gas  density  ratio  is  calculated  using  the  density  of  the  gas  in 
front  of  the  shock  wave  (the  initial  driven  section  density  pj).  Assuming  the  liquid 
density  is  constant  the  only  variance  in  the  ratio  is  due  to  range  of  pj. 

Liquid  /  Chamber  Density  Ratio  =  — 

P; 

Let  T  =  liq/gas  density  ratio  — .  Assuming  liquid  density  is  constant  for  every  test 

Py 

condition,  then: 

'^low  ~  Pl^Pglhigh  ^igh  ~  Pl^Pgllow 

9.  In  this  LRE  the  dynamic  pressure  ratio  ranges  from  0.07  to  6.64.  The  equation  for 
liquid/gas  dynamic  pressure  ratio  is  given  as:  [p^^lpg^y Knowing  the 
liquid  density,  gas  properties,  and  the  range  of  dynamic  pressure  ratio  values,  we 
calculate  a  range  of  low  to  high  jet  velocity,  Vp  that  would  result  in  matching  the 
liquid  injection  condition  to  simulate  LRE  conditions.  It  is  found  that  the  lowest 
value  of  jet  velocity  occurs  when  the  0.0374  is  used  as  the  dynamic  pressure  ratio 
and  Pg2iow  is  used  as  the  gas  density.  Similarly,  the  highest  value  of  Vj  is 
calculated  using  the  extreme  value  of  the  above  parameters. 

~  PUq  ^liq 

P«2  ''n 

therefore, 

^'9 -lew  ^low 

10.  Calculate  the  viscosity  ratio:  [vi/Vg2].  Since,  v  =  |i/  p,  the  viscosity  ratio  can  be 
expressed  as  Vi/Vg2  =  [pipg2/PgPl]-  Therefore,  due  to  dependence  on  the  range  of 
pg2,  there  will  also  be  a  low  to  high  range  of  liquid/gas  viscosity  ratio. 

11.  The  liquid  Reynolds  Number  (Re  =  piVi/pi)  is  also  calculated.  Again  due  to  a  low 
and  high  range  of  liquid  velocity,  a  range  of  Reynolds  number  would  have  to  be 
calculated. 

Examples  of  a  test  matrix  spreadsheet  are  shown  in  Table  A.1  -  A.4. 
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Table  A.l.  Test  Matrix  Spreadsheet  for  Methanol 
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Table  A.2.  Methanol  Test  Matrix  Spreadsheet  for  Various  Jet  Diameters 


Table  A.3.  Test  Matrix  Spreadsheet  for  Freon-22 


I 
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Table  A.4.  Freon-22  Test  Matrix  Spreadsheet  for  Various  Diameters 
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Appendix  B:  Jet  Displacement  Plots 


T 


Test  Wave  Type  Strength  (%)  Duration  (ms)  Diameter  (nun) 

1  N  34  5.5  0.8 

Figure  B.  1 .  Dimensionless  downstream  displacement  of  a  liquid  jet  for  Test  1 


Test  Wave  Type  Strength  (%)  Duration  (ms)  Diameter  (mm) 

2  N  57  5.5  1.0 

Figure  B.2.  Dimensionless  downstream  displacement  of  a  liquid  jet  for  Test  2. 
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Test  Wave  Type  Strength  (%)  Duration  (ms)  Diameter  (mm) 

3  N  38  2.3  1.04 

Figure  B.3.  Dimensionless  downstream  displacement  of  a  liquid  jet  for  Test  3. 


0  0.5  1  1.5  2  2.5  3  3.5 

T 

Test  Wave  Type  Strength  (%)  Duration  (ms)  Diameter  (mm) 

4  N  38  2.3  1.41 

Figure  B.4.  Dimensionless  downstream  displacement  of  a  liquid  jet  for  Test  4. 
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Test  Wave  Type  Strength  (%)  Duration  (ms)  Diameter  (mm) 

5  N  62  2.3  1.41 


Figure  B.5.  Dimensionless  downstream  displacement  of  a  liquid  jet  for  Test  5. 


Test  Wave  Type  Strength  (%)  Duration  (ms)  Diameter  (mm) 

6  S  38  infinite  1.41 

Figure  B.6.  Dimensionless  downstream  displacement  of  a  liquid  jet  for  Test  6. 
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Test  Wave  Type  Strength  (%)  Duration  (ms)  Diameter  (mm) 

7  S  65  infinite  1.41 


Figure  B.7.  Dimensionless  downstream  displacement  of  a  liquid  jet  for  Test  7. 


Test  Wave  Type  Strength  (%)  Duration  (ms)  Diameter  (mm) 

8  S  38  infinite  1.04 

Figure  B.8.  Dimensionless  downstream  displacement  of  a  liquid  jet  for  Test  8. 


T 


Test  Wave  Type  Strength  (%)  Duration  (ms)  Diameter  (mm) 

9  S  254  infinite  1.04 

Figure  B.9.  Dimensionless  downstream  displacement  of  a  liquid  jet  for  Test  9. 
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Appendix  C:  Weber  -  Reynolds  Number  Time  History  Plot 

The  Weber  -  Rejmolds  number  time  history  for  each  test  has  been  plotted  in 
Figures  C.l  -  C.9.  The  following  equations  were  used  in  determining  these  values: 


Weber  Number: 


Reynolds  Number: 


We  = 


<J 


Re  = 


p  V  D. 

rs  S  J 


The  gas  velocity,  Vg  .and  density,  pg ,  were  obtained  from  the  pressure  time  history 
plot  and  the  standard  shock  tube  relations. 

The  Weber  Number  characterizes  the  surface-tension  effect  It  is  the  ratio  of  the 
dynamic  force  of  the  gas  to  the  surface-tension  of  the  liquid.  The  Reynolds  number  is  the 
ratio  of  the  inertial  force  to  the  viscous  forces  of  the  gas. 
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Gaseous  Weber  and  Reynolds  Number  Time  History  for  Test  1 


Figure  C.2.  Gaseous  Weber  and  Reynolds  Number  Time  History  for  Test  2 


jequuriN  jeqeM 


Figure  C.3.  Weber  and  Reynold  Number  Time  History  for  Test  3 


e  C.5.  Weber  and  Reynold  Number  Time  History  for  Test  5 
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Figure  C.6.  Weber  and  Reynold  Number  Time  History  for  Test  6 


Reynold  Number  Time  History  for  Test  7 


History  for  Test  8 


Number  Time  History  for  Test  9 


